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SEMIGEOSTRPHIC EQUATIONS IN PHYSICAL SPACE WITH 
FREE UPPER BOUNDARY 

JINGRUI CHENG 


Abstract. We define various notions of Lagrangian solution in physical space 
for 3-d incompressible semigeostrophic system with free upper boundary under 
different conditions for initial data, then prove their existence via the mini¬ 
mization with respect to a geostrophic functional, generalizing the the work of 
go and Eli to the case of free upper boundary. As a byproduct of our proof, 
we obtain the existence of measure-valued dual space solutions when the initial 
measure uq € ^(R- 3 ) and is supported on {—j < ys < —<5} 


Contents 

1. Introduction 

2. The study of the functional E v (h, 7) 

2.1. The case when v has bounded support 

2.2. Generalization to v with unbounded support 

3. Existence of Lagrangian solutions 

3.1. Basic definitions and main result 

3.2. Lagrangian flow in dual space 

3.3. Lagrangian flow in physical space 

4. Existence of relaxed Lagrangian solution with general initial 

data(VPo G L 2 ) 

4.1. Definition of generalized data and main result 

4.2. Measure-valued solution in dual space 

4.3. Existence of relaxed Lagrangian solution with generalized data 

5. Appendix 
Acknowledgments 
References 


1 

5 

J£ 

17 

24 

24 

3C 

32 


34 

34 

36 

39 
43 

40 
46 


1. Introduction 

The Semi-Geostrophic system (abbreviated as SG in the following)models large- 
scale atnrospheric-ocean flows, where large scale means the flow is rotation-dominated, 
J.-M. Benamou and Y.Brenier [4: proved the existence of solutions to incompress¬ 
ible SG system in a fixed domain D C R 3 in the so-called dual space formulation, 
which is a formal change of variable. Under the dual space formulation, the SG 
system can be written as a transport equation coupled with Monge-Ampere equa¬ 
tion. Mike Cullen and W.Gangbo [6] considered the free boundary case in 3-D, but 
with additional assumption that the potential temperature is constant. Under this 
additional assumption, the system can be rewritten as a system in 2-D, the so-called 
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Semi-geostriophic Shallow Water system. The authors then proved the existence 
of dual-space solutions with initial dual density in L p (p > 1). The existence of 
solutions in the original physical variables is first proved by Cullen and Feldman 
[5] , in the Lagrangian formulation of the physical system, for both fixed boundary 
SG system and shallow water system with the same assumption on the dual den¬ 
sity as above, and this assumption amounts to some strict convexity condition for 
the modified pressure. Then in mm, Feldman and Tudorascu put forward a 
more general notion of physical solutions, to allow for more general initial data, in 
particular, they proved the existence of general measure valued solutions in dual 
space. 

In this work, we consider the incompressible SG system in a 3-D domain with free 
upper boundary, but without the constancy assumption made in [ 6 ]. The existence 
of dual space solutions has been proved in [7] by invoking the general theory of 
Hamiltonian ODE established in [2]. Here we prove the existence of Lagrangian 
solutions in physical space, generalizing the work of [5] and m The main difficulty 
involved is the more complicated geostrophic energy in our situation since it involves 
the unknown free boundary profile. The 3-D SG system with free upper boundary 
has the following form: 

(1.1) D t ( uf,u|) + ( u 2 , ui) + (d Xl p,d X2 p) = 0 

(1.2) V - u = 0 

(1.3) D t p= 0 

(1.4) S7p= (uf,-uf,-p) 

(1.1)-(1.4) hold in Q hl where D t = d t + u ■ V, and Sl/j = {(xi, x 2 , £ 3 ) £ fl 2 x 
[ 0 ,oo )|0 < x 3 < h(t,xi,x 2 )} 

Here fi 2 C R 2 is a bounded convex region. h(t,x i,x 2 ) > 0 describes the un¬ 
known free upper boundary. In the above p is the pressure, u is the velocity, and 
p is the density. 

Of course we need to prescribe suitable boundary and free boundary conditions. 
We require that no flow can penetrate the fixed boundary, the pressure at the top 
is a constant which without loss of generality we take to be zero. This is given 
by (1.5) and (1.6) below respectively. (1.7) means that no flow penetrates the free 
boundary, instead the fluid particles move with the flow. 

(1.5) u • n = 0 on dflh — {£3 = h} 

( 1 . 6 ) p(x ll x 2 ,x 3 ) = 0 on {x 3 = h} 

(1.7) d t h + n x d xi h + u 2 d X2 h = u 3 on {x 3 = h} 

We remark that (1.5) and (1.7) combined express conservation of mass , and are 
formally equivalent to 

d t o- h + V • (ucr h ) = 0 

where 


(1.8) 


cth{x i,x 2 ,x 3 ) = X'n h (xi,x 2 ,x 3 ) 
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Now we put 
(1.9) 

then the above system can be written as 


P(t, x) = p(t , x) + ~{x\ + x\) 


( 1 . 10 ) 

( 1 . 11 ) 


A (VP) = J(VP — x) 
V • u = 0 


In the above 

/ 0 -1 0 

J= 1 0 0 

\ 0 0 0 

The free boundary condition for P is 

1 


( 1 . 12 ) 


P(x i,X2,x 3 ) = ^(xi +xl) on {;c 3 = h} 


The geostrophic energy is 

(1.13) 

E = [ i((uf ) 2 + {u 9 2 ) 2 ) + px 3 dx = f ]-[(dxP - xx) 2 + {d 2 P - x 2 ) 2 } - d 3 Px 3 dx 
z Jn h z 

By Cullen‘s stability principle, the function P above should be convex, and VP 
should minimize above functional among all possible rearrangement of particles. To 
make this precise, we are motivated to consider the following minimization problem. 

(1.14) E v (h, T)= f ^[(xx-Txix)) 2 + (x 2 -T 2 (x)) 2 ]-x 3 T 3 (x)dx 

J Ooo 1 

Here v £ P(R 3 ) is fixed, T^oy,, = v and we require that the actual free boundary 
profile h and the Borel map VP minimizes E„(h, T) among all pairs of (h, T) such 
that h > 0, continuous, h = 1 , and Tjay = v. 

Also we recall that d X3 P = —p, where p represents the density, so it‘s reasonable 
to assume the convex potential P should satisfy — ^ < d X3 P < —5, or we require 
supp v C R 2 x [—4, —5], More generally, we can consider 

(1.15) E v {h, 7 ) = J [i[(xi - yi) 2 + (x 2 - y 2 ) 2 ] - x 3 y 3 }d^(x,y) 


Here 


h > 0, / hdx\dx 2 = 1 7 € r(ay,u) 

J 0,2 


Here T(ay,u) denotes the set of measures on the product space Hoc x R 3 with 
marginals ay and v. Here floo = fi 2 x [0,oo). This is the relaxed Kantorovitch 
problem of (1.14). For a fixed ay, if 7 £ T(ay,u) minimize (1.15), then 7 is called 
an optimal transport plan of (1.15). If the optimal plan has the form 7 = (idx T)jj<jy 
for some Borel map T, then this problem reduces to (1.14). 

In the following, we will denote the cost function: 

(1.16) c(x,y) = ^[(aq - yx) 2 + (x 2 - y 2 ) 2 ] - x 3 y 3 = ^(xf + x\ + y\ + y\) - x ■ y 

Here x = (xx,x 2 ,x 3 ) e Hoc and y = ( 1 / 1 , y 2 , y 3 ) € R 3 . 
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Define v = VPjcr^, where P is the function appearing in (1.10), then formally 
it satisfies the following continuity equation. The derivation can be found in [7|, 
section 2. 

(1.17) d t v + V- (i/w) = 0 in [0,T] x A 

(1-18) w (t,y) = J(y-\7P*(t,y)) m[0,T]xA 

(1.19) (h, VP) minimizes E v (h,T) among all pairs such that Tjcr^ = v 

Pelloni et al proved in [7J the existence of weak solutions to above SG system in 
dual variables, (1.17)-(1.19), using the theory of Hamilton ODEs developed in [1| . 
In this paper, we will reprove this result using a more straightforward time-stepping 
argument. Then we prove the existence of solutions to physical equations in La- 
grangian variables when the dual density v defined above is absolutely continuous. 
Finally we define a notion of relaxed Lagrangian solution similar to m and prove 
its existence. This allows us to deal with the case where v £ P 2 (R 3 ) (probability 
measures with finite second order moment). 

Here are some notations and terminology which will be used thoughout this 
work. In the following, Hoc represents the region H 2 x [0, 00), and fin = ^2 x [0, H], 
Given h : H2 —>• R+, we denote flh = {x £ Hoo |0 < x% < h{x 1,0:2)}, and 0h(x) = 
Xn h {x). We identify an absolutely continuous measure (with respect to C 3 ) with 
its densities. 

Suppose A,Bc R 3 , and we have two functions f(x),g(y) defined on A and B 
respectively. We say /, g are convex conjugate to each other over A and B if the 
following holds. 

fix) = sup ix ■ y - g(y)) x £ A 

y£B 

and 

g(y) = sup (a: • y - /( x)) y £ B 

xeA 

In the sequal, let p > 1, we denote P p (R 3 ) to be the set of probability measures 
in R 3 such that f R3 \y\ p dy < 00. For v £ P p (R 3 ), denote M p {v) = f \y\ p 

In the following, V 2 mans the gradient with respect to the first 2 variables only, 
namely V 2 U = (d\v,d 2 v) 

Given p, v £ P p (R 3 )(p > 1), we define the p-Wasserstein metric to be 
W p {p,v)= inf f \x-y\ p d~/{x,y) 

7er( M ,y) y R 3 xR 3 

Here T{p,v) denotes the set of probability measures on R 3 x R 3 with marginals 
p. v respectively. W p (p, u) is indeed a metric on P p (R 3 ), see [3] chapter 7. 

To conclude this section, we briefly describe the plan of this paper. 

In section 2, we study the geostrophic functional and it‘s dual problem in detail, 
and establish various properties of the optimizers which will be used later on. Then 
in the case when the dual density v £ L q , for some q > 1 with compact support, we 
follow [5] to establish the existence of weak Lagrangian solutions, using the theory 
of Lagrangian flows generated by BV vector fields developed in [T], see Theorem 3.9 
in section 3.3. In the case when v is singular and may have unbounded support, we 
generalize the notion of weak Lagrangian solutions and prove their existence with 
suitable initial data, see Theorem 4.6. As a byproduct, we obtain the existence 
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of measure-valued dual space solutions when the initial dual density vq 6 ^(R 3 ) 
with support contained in R 2 x [—4, — <5] for some S > 0, see Corollary 4.13. 


2. The study of the functional E v (h,"i) 


2.1. The case when v has bounded support. In this section, we study the 
functional involved in the geostrophic energy and the associated dual problem, 
prove basic properties such as unique existence of optimizers. Finally we give an 
alternative proof of dual space existence result using time stepping since later on 
we will need some regularity properties of dual space solutions which are not so 
clear in Hamiltonian ODE approach as was done in [TJ. 

We study the property of the functional 

(2.1) E„(h, 7 ) = / c(x, y)d r ){x, y) = f [\{x\ + x\ + y\ + y\) - x ■ y\d^ 

JQ ooXA JQ ooXA ^ 

where (h, 7 ) £ M. v v £ ^(R 3 ) and supp v C A C R 2 x [— |, — d] and compact. Here 
let‘s choose A = Bb{ 0) x [— —5] and assume for technical reasons C Bjj( 0) 

where 


Aty = {(h,"f)\h > 0, continuous, / hdx\dx 2 = 1,7 G r(cxfc,,u)} 

Jn 2 

and also the functional 

( 2 . 2 ) 

J?{P,R)= [ [\(yl + yl)-R{y)]v{y)dy+ inf / {hxj+x%) - P(x)\a h {x)dx 
J A 2 0 <h<H J Qoo 2 

where 

P{x) + R(y ) > x ■ y Vx G x [0, H] \/y G A 
Call the collection of all such pairs satisfying the above condition to be M v . 

We will also consider the untruncated version, namely 

(2.3) J V (P,R) = [ \^{yl+yl)-R(y)\v{y)dy+mt [ [\{x\+xl)-P{x)}a h {x)dx 
J A 2 h >°J 2 

where we require 

R G L 3 (du) P G L 1 {VLk) VK > 0, and h G L 1 (H 2 ) Pah G L 1 (H 00 ) 


and 

P{x) + R(y) > x • y Vx G floo y G A 

We will see later that J,y(P,R) is dual to E v (h^) in the next subsection. Also the 
study of the dual problem J V (P, R) will help with proving uniqueness of minimizers 
of E v (h, 7 ), since the geostrophic functiona^-E^ (h, 7 ) in our case does not seem 
to have strict convexity as in mm- This was first noted in [7] and the idea of 
considering a dual variational problem is inspired by [5]. 

In the case when v has bounded support, A can be taken to be bounded, and 
later on(in this subsection) we will show the h which assumes the infimum in (2.3) 
has a universal bound in L°° (depending only on the data of the problem) and so it 
will be equivalent to solving the truncated problem J^(P, R) if one takes H large 
enough depending only on 5, H 2 and A. 

Suppose d X3 P(x) < — 5 , let‘s define 

f S 1 

n P (xi,x 2 ,s) := / [~(xl +X 2 ) - P(xi,x 2 ,x 3 )]dx 3 
Jo ^ 
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It‘s easy to see for fixed (£ 1 , 0 : 2 ) £ O 2 , the function s 1 —> n p (£i,£ 2 ,s) is uni¬ 
formly convex and so there exists a unique s* where lip achieves minimum on 
[0,oo). We define this function to be hp(x i,£ 2 ). We also define hp{x i,£ 2 ) to be 
the unique s* £ [0,//] where lip achieves minimum on [0,1/]. Notice by convex¬ 
ity,one has hp = min (hp, H). 

Remark 2.4. Whenever hp(x i,£ 2 ) > 0, we must have 

P(x 1 ,x 2 ,h P (x 1 ,x 2 )) = ^(£1 +£ 2 ) 

Otherwise 

P(£!,£ 2 ,0) < -(£? + £ 2 ) 

Conversely, if P is defined on Q 00 and h satisfies above condition, then we will also 
have h = hp 

Remark 2.5. IPs easy to see in the situation of J„(P,R) 

(2.6) inf [ [~{x\ + x%) - P(x)}ah{x)dx = f [~{x\ +x%) — P(x)]ah p {x)dx 
and in the situation of J„(P, R) 

( 2 - 7 ) nJ&v f \k{x\ + xl)-P{x)]a h {x)dx= j []-(xl+xl)-P(x)\a h H(x)dx 
0<h<Hj nH 2 Jn H 2 p 

Now we prove the following 

Lemma 2.8. Suppose there exists a sequence d X3 P n , d X3 P < —5, and P n —>• P 
uniformly on tt 2 x [0,i/] for each H > 0, then we have hp n (xi,x 2 ) —> hp(xi,x 2 ) 
uniformly on Vt 2 . If P n ,P satisfy the same condition but is only defined on fl 2 x 
[0,i/o], then hp° —> hp° uniformly on Vl 2 

Proof. Denote h n = hp°, h = hp° for convenience. First we show that {h n } is 
uniformly bounded. Indeed, 

t)(£? +£ 2 ) = P n {x 1 ,x 2 ,h Pn (x 1 ,x 2 )) < —Sh n (xi,x 2 ) + P n (xi,x 2 ,0) 

Put e n = sup xen2X[0j p] |P„(£) - P(£)|, then we have 
If hp(xi,x 2 ) = 0, then 

i(£’l+£ 2 ) > P(x 1 , £ 2 , 0 ) > -e„ + (P„(£ 1 ,£2,0)-P n (£ 1 ,£ 2 ,/l„(£ 1 ,£ 2 )) + ^(£? + £ 2 ) 

> —e„ + Sh n (x i,£ 2 ) + i(£? + £ 2 ) 

If hp(xi,x 2 ) > 0,and h n > h then 

0 = P n (xi,X 2 ,h n (x i,£ 2 )) - P(£i,£ 2 ,/l(£l,£ 2 )) 

= P n (xi,x 2 ,h n ) - P(xi,x 2 ,h n ) + P(xi,x 2 ,h n ) - P(xi,x 2 ,h) < e n - 5{h n - h) 
The other case h n < h can be dealt with similarly. □ 

We also need another lemma which gives control over the absolute bound for the 
maximizing sequence. 
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Lemma 2.9. Suppose there exists constant K > 0 and P : Qh —> R with P(xi,, x\, 0) 
0 where £ fl 2 and P,R are convex conjugate over the domain fin and A 

respectively, such that for some A > 0 

—K < jf(P-A,P + A) 

and 

|V 2 P| < K - i < d x ,P < -6 

then 

|A| < Ci(K, A, Cl 2 , H) 

as long as H > 

Proof. First from our assumption on P , we obtain 

\P(x)\ < Kdiam fi 2 H—- ; = C_ 1 x £ f Iff 
0 

Since (P, R) are convex conjugate, one has 

R(y) = sup (x ■ y - P(x)) > x\yi + x\y 2 > -max|x| -max(|j/i| + \y 2 \) := -Co 
xen H J/eA 


From the definition of J H , by taking h = 0, one has 

-K < f [hyl + yl) - R(y) - A \v(y)dy < max(|j/i| + \y 2 \) 2 + C 0 - A 
A 2 v^ A 

So 

A < C 0 + K + max(|j/i| + \y 2 \) 2 

ye A 

On the other hand we take h = we then have 

~K < f \hy 2 + yl) - R(y) - A \v(y)dy + [ + a|) - P(x) + \\dx 

J a 2 ^n 2 x[o I? %] 2 


< maxflj/il 
yeA 


+ |2/2|)“ + C'_i + C'o + 


max \x\ 

o 2 


+ A 


□ 


Now we prove the existence of a pair of maximizer of Jff, using a standard 
compactness argument. 

Theorem 2.10. Suppose H > ^i 2 ^ > then the variational problem J^f has a 
maximizer (. P,R ), where P and R are convex conjugate to each other over Qh and 
A 


Proof. We choose a maximizing sequence (. P n ,R n ), without loss of generality, we 
can assume they are convex conjugate over the domain Qh and A, then their 
derivatives are uniformly bounded(with a bound depending on fl 2 ,A , H) and P 
satisfies A- < d X3 P < —6. 

Now the functions 

P n :=P n ~P n (xl,xl 0) R n ~ R n + P n (a$ , X \ , 0) 

satisfy the assumption of previous lemma with A = —P n (x\,x 2 ,0), we can then 
conclude 


|P n (*i,a:5,0)| <C 
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Since the derivatives are bounded independent of n, we get P„ is uniformly 
bounded independent of n. Also since P n ,R n are convex conjugate, we see R n are 
also bounded uniformly independent of n. 

Now we can apply Arzela-Ascoli to get a pair (P, R) which is also convex conju¬ 
gate over Hp and A. 

Also by lemma 2.8, we also have 


h p n h P 

uniformly on Sl 2 . Since we also have 

P n —> P R n —> R uniformly 

One can see 

J?(P n ,R n )^J?(P,R) 

Hence (P, R) is a maximizer. □ 

Lemma 2.11. Suppose for some H > 0, the variational problem Jff has a pair of 
maximizer (P, R), such that P(x) = sup yeA (:r • y — R{y)) and put h = hp , then 

\7P$cr h = v 

in particular, 

[ h = 1 
Jn 2 

Proof. Let g(y) £ Cb(A),5 > 0, define 

Rs(y) := R(y) + 5g(y) Ps{x) := sup(z • y - Rs(y)) 

yeA 

Since (P, R) is a pair of maximizers, we have 

J H (P s ,Rs ) < J h (P,R) 

i.e, we have 

-5 [ g(y)du(y) < inf [ [\(xl+xl)-P{x)]dx- inf f \Ux\+xl)-P s {x)}dx 
J A 0 <h<H J Qh 2 0 <h<H J Qh 2 


< / dxidx 2 / [P s (xi,x 2 ,x 3 ) - P(x 1 ,x 2 ,x 3 )]dx 3 

Jn 2 do 


Here 


hs(x i,x 2 ) = hp 5 (x i,x 2 ) 

Note that np^ (xi, x 2 , s) achieves minimum over [0,P] at hs■ Therefore hs 
achieves the second infimum above. 

Then we notice that since Rs —> R uniformly on A, we have Ps —> P uniformly 
on Qh- Hence by Lemma 2.8, we have 

hs(xi,x 2 ) h(x\,x 2 ) uniformly 

Suppose P is differentiable at x, and let ys € A be the point such that 

P s (x) =x-y s - R s (y s ) 

then we have 

ys —> VP(a;) as S —> 0 
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Also notice that 


~g(VP(x)) < < -g{y s ) 


By letting S —» 0,we obtain 


- [ g{y)v{y)dy < - [ g(VP(x))dx 


J A J n h 

Replacing g by — g , we are done. 

Now we prove the function hp obtained above is Lipschitz. 


□ 


Proposition 2.12. Let (P, R) be a pair of convex conjugate maximizers over the 
domain fi 2 x [0, if] and A respectively, then hp (x± , 2 : 2 ) is Lipschitz, with a Lipschitz 
constant depending only on 6, A, fi 2 (not H ) 

Proof. Let's denote hp(x 1 , 2 : 2 ) = h. The proof is based on the following facts. 

(■ i)P(xi,X 2 ,h(xi,X 2 )) = 5 ( 2:1 + x\) whenever 0 < h < H 
( ii)d X3 P < —5, and |V 2 P| < C. Here C depend only on A 
Pick ( 2 : 1 , 2 ; 2 ), ( 21 , Z 2 ) G n 2 , without loss of generality, we can assume h{z\, 2 2 ) < 
h(x 1 , 2 : 2 ) 

If 0 < h(zi, Z 2 ) < h(x 1 , 2 : 2 ) < if, then we have 

P(xi,X2, h{x\,X2)) = ^(xf + a|) P(z 1 ,z 2 ,h(zi,z 2 )) = 7j(z? + z%) 

Thus 


-C{\x 1 -zi\+\x2-z 2 \) < P(xi,x 2 , h(xi,x 2 ))-P(x 1 ,X2,h(zi, Z 2 ))+P(xi,x 2 , h(z 1, 2 2 ))+f , ( 2 i, 22, h(zi, z 2 )) 
< -S(h(xi,x 2 ) - h(zi,z 2 )) + C(\xi - 2i| + |x 2 - 2 2 |) 

Next consider 0 < h(2i,2 2 ) < /i(x i, 2 : 2 ) = if, then we have 

P(xi,X 2,H) > ^{xl+xl) P(z 1 ,z 2 ,h(z 1 ,z 2 )) = ^(21 +2I) 

Thus 


— C'daii - 2 i| + |x 2 - 2 2 |) < P(x 1 ,x 2 ,H) - P{z 1 ,z 2 ,h{z 1 ,z 2 )) 

< C(\xi - 2 i| + | 2 J 2 - 2 2 |) - S(H - h{z\ , 2 2 )) 

Notice that if h( 2 i, 2 2 ) = 0, then P(zi,Z 2 , 0 ) < \(z\ + z$), see Remark 2.4, so 
other cases can be dealt with in a similar way. □ 

Corollary 2.13. Suppose if > (such that jff{P,R) has maximizers by The¬ 

orem 2.10), and let ( P,R) be convex conjugate maximizers defined on Qh and A 
respectively, then there exists a constant C = C(diam fl 2 ,<5, A), such that 

hp(xi,x 2 ) < C 


Proof. By Lemma 2.11, we have 



so there exists ( 21 , 22 ) G fl 2 , such that /i( 2 i, 2 2 ) < 

By above corollary, we know h is Lipschitz, and since 

P{x 1 , 2 ? 2 , h) = — ( 2 ;^ + a: 2 ) whenever h > 0 
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SO 


Thus 


So 

(2.14) 


Vh(xi, x 2 ) 


-——{x 1 -d 1 P,x 2 

U X3 1 


d 2 P) 


|V/i| < 


maxo 2 \x\ +max A (|(/i| + \y 2 \) 

8 


\h\< 


c 2 (n 2 ) 


+ 


2 maxn 2 \x\ +max A (|j/i| + I 2 / 2 I) 


• diam fl 2 


□ 


As an easy consequence, we deduce: 

Corollary 2.15. Suppose H > Co, here the Co is the constant on the right hand 
side of (2. If), if{P,R) is a pair of maximizer of jff(P,R), convex conjugate over 
n H and A, then 0 < hp < H. 

Corollary 2.16. Suppose H satisfy the same condition as in previous corollary, 
then there exists another constant C = C(S, fl 2l A, H), such that if (P, R) is a convex 
conjugate maximizer of jff(P,R), we have |P(x)| < C, V x € £Ih 

Proof. Fix (xi,x 2 ) £ Sl 2 , if h(xi,x 2 ) > 0, we then have 

7 }( x i + x l) < P(xi,x 2 ,0) 

On the other hand 

1 f h 

-(xl +xl) = P(x 1 ,x 2 ,h(xi,x 2 )) = P(xi, x 2 ,0) + / d X3 P(xi,x 2 ,s)ds 

* Jo 

> -\h{x 1 ,X 2 ) + P(x±,x 2 , 0) 

0 

Thus by Corollary 2.13, we know 

\P(x±, x 2 , 0)| < C such that h{xi,x 2 ) > 0 
where C has the said dependence. 

Now notice that since by assumption (P, R) are convex conjugate over the domain 
Qh and A respectively, we know 

dP(n H ) c A 

Hence 

|VP(*)| < C{A,S) a.e 

Fix ( 2 : 1 , 22 ) such that h(z\,z 2 ) > 0, then for x £ fin, we have 

P(x 1 ,X 2 ,X 3 ) = P(zi,Z2,0)+ / VP(t(zi, Z2,0) + (l-t)(xi,X 2 ,Xs)y(xi—Zi,X 2 -Z 2 ,X 3 )dt 

Jo 

The result follows easily. □ 


In the sequel, we will always assume 


(2.17) 


H > 


£ 2 (^ 2 ) 


2 maxfi 2 \x\ +max A (|yi| + \y 2 \) 

--- • dlCLTYl H2 


unless otherwise stated. 
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Theorem 2.18. Suppose H is as in (2.11), v £ P(R 3 ), with supp v C A then we 
have 

(' i)E v (h , 7 ) > (P,R), for any (h, 7 ) £ M. v and (P, P) satisfying P(x)+R(y) > 

x ■ y \/x £ fin Vy £ A 

(ii)Suppose (P, P) is a convex conjugate maximizer of J H (P, R), then we have 
equality in (i) if and only if 7 = (id x VP)jot, and h = hp 
(in)In the situation of (ii), and if v « C?, we also have 

VPju = (jfi 

Proof. First we prove (i). 

Without loss of generality, we can assume (P, R) be a pair of convex conjugate 
maximizers of jff (P, R), and we can naturally extend P to be defined on such 
that (P, R) are again convex conjugate. Indeed, 

P(x) = sup(a; • y - R(y)) 

v&A 

and we define P(x) for x £ floo by the same formula. Then one can check 
R(y) = sup (x ■ y — P(x)) = sup (x ■ y — P(x)) 

x^Cloo 

This implies in particular 

P(x) + R(y) > x ■ y x £ y £ A 
Therefore, assuming (h, 7 ) £ JA, we can write 


c(x,y)d'y(x,y) > 


[\(yl+yl)-R(y)\v(y)dy- 


\^(x\ + x%) - P(x)\a h (x)dx 


> 


- } (yl + yl)-R(y)]iz(y)dy + \^(x\ + x%) -P(x)\a h n(x)dx = J H (P,R) (*) 

Here only the second ineciuality requires some explanation. 

By Corollary 2.15 and by strict convexity of np(xi, X 2 , s) in s, we see that 
s 1 —»• TVp(x\, X 2 , s) attains minimum over [0,oo) at hp(x 1 ,^ 2 ). So we have 

[ \t.(x\ + x\) - P(x)\ah(x)dx = I Hp(xi,x 2 ,h(xi,x 2 ))dxidx 2 

Jn a o 2 " Jn 2 

> f U p (xi,x 2 ,hp(xi,x 2 ))dxidx 2 = f [-(x\+x\) — P(x)]a h H(x)dx 

Jn 2 Jn^ 2 p 

It‘s easy to see above inequality takes equality if and only if h(x 1 , 212 ) = hp (x\,X 2 ). 
Up to now,we proved (?) 

Now we can prove (ii). Suppose we have P y (/i, 7 ) = j(f(P,R ), then both in¬ 
equality in (*) must be equality. Thus we know h = hp and we know from Corollary 
2.13 that supp ah C Up. The first inequality in (*) takes equality if and only if 
P(x) + R(y) = x ■ y 7 — a.e (x, y), or equivalently y £ dP(x) 7 — a.e (x, y). So we 
obtain 7 = (id x VPjjcr/j. Conversely, if h = hp and 7 = (id x VPjjd/,, then by 
Lemma 2.11, (h, 7 ) £ M„, also we have E v (h , 7 ) = jff(P,R). 

To see (Hi), note that by (ii), we have 

VPflcr^ = v 

And since v « C?, we know 

VP o VP(x) = x Oh a.e 
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Indeed put 


E = {y e A|VP is not defined at y} 


then 

v{E)=a h ((yP)-\E)) = 0 

And if VP is defined at x, VP is defined at VP(x), then one has VP(VP(x)) = x. 
So 


O'/I = VPttVPjcrft = VP#!/ 


□ 


Now we can prove the unique existence of the variational problem E v (h, 7 ). 

Corollary 2.19. Suppose H is as (2.17), then 

(i) E v (h , 7 ) ftas a unique minimizer (h, 7 ); 

(ii) Suppose (Pq, Pi),(Pl, Pi) are maximizers of Jff(P,R), convex conjugate over 

Qh and A, then h.p Q = := h, and Pq = Pi on 12/, 

Proof. Existence of at least one maximizer of (P, P) has been proved in Theorem 
2.10. Theorem 2.18 (ii) gives the existence of at least one minimizer of E„(h, 7 ). 
Now we show uniqueness. 

First we fix a maximizer of jff(P,R), say (Po,Po), which is convex conjugate 
over 1 Ip and A. If (ho, 70 ),(/ii, 71 ) are both minimizers of E v (h, 7 ), then we have 
by previous theorem 

ho = hi = hp 0 (xi,x 2 ) := h 0 , and ji = 70 = (id x VPojjtT/j 
This proves the uniqueness of the minimizer of E v (h, 7 ). 

To see the uniqueness of J H (P, P) on 12 h 0 , say both (P 0 , Po),(Pi, Pi) are maxi¬ 
mizers. Let /io = and 70 = (id x VPojjcr/jo. Then we know by Proposition 2.12 
that ho is Lipschitz, also we know by Theorem 2.18(iz) that h = hp = hp Q 

Also by uniqueness of minimizer for E v (h, 7 ) already proved above, we know 
from above theorem (ii) 

(id x VPojjCTfto = (id x VPijjjcr^ 

This implies 

VP 0 = VPi a ho a.e 

Let U C 122 to be a connected component of the open set {(x\,x 2 ) £ fl 2 \ho(xi,x 2 ) > 
0}, then the set Uh ■= {(x\, x 2 , X 3 )\(x\, x 2 ) £ U, 0 < X 3 < ho(xi,x 2 )} is connected 
and so we have 

Pq = Pi + Cu on Uh 

But we also have 

Po(xi,x 2 ,h 0 (xi,x 2 )) = Pi(xi,x 2 ,h 0 (xi,x 2 )) = ^(xf + x\) 
for (x\,x 2 ) £ U, hence Cu = 0. So we have 


P 0 = Pi cr ho a.e 


□ 
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Corollary 2.20. Let A = Bd{ 0) x [—^,—5] and assume also that LI? C Br(0), 
H is chosen as (2.17). Then there exists a unique maximizer (P 2 ,P 2 ) of jff(P,R) 
with the following properties: 

(i) (P 2 ,R 2 ) are convex conjugate over both = 0}, A and Q,fj, A. 

(ii) P 2 (xi,X 2 , 0) = \(x\ + x 2 ) whenever ho(xi,x 2 ) = 0 

Proof. The uniqueness is easy to see, since by (i), we must have 
(2.21) R 2 {y)= sup (x-y-P 2 (x)) My G A 

X&O.h 0 U{ a; 3=0} 

The previous corollary shows that for any maximizer (P, R) convex conjugate over 
LIh and A, the value of P must agree on Oft, while point (ii) takes care of the part 
where X 3 = 0 and h(xi,x 2 ) = 0. Hence R 2 is uniquely defined by (2.21). 

And 

P 2 (x) = sup (31 • y - P 2 (y)) Vcc G Q. h 
ye A 

because (P 2 ,R 2 ) are assumed to be convex conjugate over f2#, A, so P 2 is uniquely 
defined over LIh- 

The existence of such a maximizer is more technical and is proved in the appen¬ 
dix. □ 


To conclude this section, we notice the following stability results. Their proofs 
are standard compactness argument, making use of the uniqueness proved in Corol¬ 
lary 2.19 and 2.20. We sketch the proofs. 

Lemma 2.22. Let v n —> v narrowly, then we have 

inf E v (h, 7) < lim inf inf E Vn {h,'y) 

(h,~t)eMv ra-s-oo (h,~i)eMv n 

Proof. Notice the uniform Lipschitz continuity and boundedness of h n proved in 
Proposition 2.12 and Corollary 2.13, and weak compactness of j n , finally make use 
of the uniqueness given by Corollary 2.19 ( i ). □ 

Lemma 2.23. Let u n —> v narrowly, then we have 

sup J I f f(F > ,R) > lim sup sup J^ n {P,R) 

(P,R)eM n->00 (P,R)^M Vri 

Proof. Choose maximizers (P„, R n ) according to Corollary 2.20, using the bounded¬ 
ness given by Corollary 2.16 and that they are uniformly Lipschitz continuous. □ 


Corollary 2.24. Suppose v n —> v narrowly, then 


lim inf E Vn (h, 7 ) 
n—too (ft,,7 )£M„ n 


inf 


K(h, 7 ) 


Proof. ” > ” follows from Lemma 2.23.” < ” follows from Lemma 2.22. Recall that 
by Theorem 2.18 (m), we have 

inf E„(h, 7 )= sup j(f(P,R) 

(P,R)£ 


□ 


Next we prove a stability result of the optimizers under narrow convergence. It 
will be useful when one proves the continuity in time of certain quantities. 
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Theorem 2.25. Suppose v n ,v £ P(A) and v n —>• v narrowly. Let (P n ,R n ), ( P,R ) 
be the unique maximizers of given by Corollary 2.20, (h n ,'y n ),(ho,'yo) be 

the minimizers of E Vn ,E u respectively, then the following holds. 

(i) h n —>• ho uniformly, 

( ii) P n —» P in W l ’ r (Q. H ), for any r < oo 

( iii) R n —»■ R in W 1,r (A), for any r < oo 

Proof. This is proved in a similar way as Lemma 2.22 and Lemma 2.23. We simply 
note that by convexity, the uniform convergence of P n ,R n implies the a.e conver¬ 
gence of VP„,VP„, and they are uniformly bounded, hence the claimed convergence 
follows. □ 

Next we apply the above obtained results to give an alternative proof for the 
existence of weak solutions in dual spaces with absolute continuous initial data, 
which was first done in [7], using Hamiltonian ODE approach. Our approach here 
is less involved. 

Now we define 

(2.26) A 0 = B Do (0) x [-^ -5] and A = B D (0) x <5] 

Based on what we proved above about the properties of the functional E„(h, 7 ) 
and jff(P,R) in Theorem 2.18, the dual space problem can be reformulated as 

(2.27) d t v + V • (y w) = 0 

(2.28) ™ = J(y-VP*(t, y )) yen 3 

(2.29) (h,(idxA7P)§(Jh) minimizes E u (h,"/) among 

In the above, P*(t,y) = sup xg Q H (a: ■ y — P(x)) for some large enough H 


Theorem 2.30. Let 1 < q < 00 , T > 0 be given and uq £ L q (Ao), where Ao,A 
defined as above, suppose also 

(2.31) 
and 


D > Dq + max |rc| • (T + 1) 
^2 


(2.32) 


H > 


+ 2m "°d I ' + 2 - D ■ diam Si, 


C 2 (n 2 ) 

Then there exists a weak solution to (2.27)-(2.29), ( h,P,R ) withu(t,-) := VP^cr^,.)) 
R = P* such that 


(■ i)t ' - > v{t, ■) £ Pac(R 3 ) is narrowly continuous with supp v C A 

v[t,-) £ L q ( A), and \\v(t, -)IU«(A) = IM')IU«(A), Vt £ [0,T] 
(ii)(P(t, •), R(t, •)) is the unique maximizer of (P,R) 
with properties (i),(ii) in Corollary 2.20 
{iii)Wi(v(ti,-),v(t 2 , •)) < C\ti - t 2 1 
(iv)P(t, •) € L°°([0,T], W 1,00 (r2 00 )) f) C([0, T], H /1 ’ r (D 00 )) 

R £ P°°([0, T], W 1,00 (A)) P) C([0, T], W 1,r (A)) Vr < 00 
(• v)Let y(t, •) = (id x VP)# 07 ,, (h, 7 ) are minimizers of E„(h,"/) and 
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h G L°°([0,T], W 1 ' oo (f) 2 ))f|C'([ 0 ,T] x n 2 ) 


Proof. Let j s : R 3 —> R be standard modifier defined by 

Mv) = ) 

First we mollify the initial data by defining 

v° 8 (y) =js*Mv) 

Then for 0 < s < 5 , we have 


2 S 

supp C B d o+s ( 0) x --) 


Next we construct approximate solutions, here we need to control the speed of 
propagation of the support of approximate solutions. Define 

D k = D 0 + ks max lari 0 < k < — 

O2 s 

Given u k := v s (ks , •) with supp n k C B D k(0) x [— |, —|], we obtain v k+1 such 
that supp v k+1 C B D k+ 1 ( 0 ) x [—|, — < 5 ] in the following way. 

Let (h k , 7 *) be the minimizer of E v k(h, 7 ), Let (P k ,R k ) the unique maximizer 
of J^ k (P, R) given by Corollary 2.20. This is possible by our assumption on H and 
D, see (2.17),(2.31) and (2.32). 

We set 7 ^ = (id x VP k )$cr h k. 

Define 

w k s (y) = J(y - VR k (y)) 
u k (y) = J(y-V(j s *R k )(y)) 

Observe that u k (y) is C°° and divergence free and we obtain v k+1 by solving 
the transport equation 

d t v s = —V • (u k v k ) in [fcs, ( k + l)s) x R 3 

u s (ks,y) = u k (y) in R 3 
and set v k+1 = v s ((k + 1 )s,y). 

In more detail, we solve the ODE 


dt 


— u s (®s (L V )) t € [ks, (k + l)s],y G R‘ 
$s(ks,y) = y 


Then 

v a {t,y) = v k (($ k )- 1 (t,y)) t G [fcs,(fc + l)s] 
since u^' is divergence free. 

Therefore 

supp v s (t, •) C -)(supp v k ) 

Now take yo G supp v k , we know 


d\**(t,yo)\ * 


k . d<S> k s 
s ' dt 


dt 


l*$l 


* k a -J($ k s -V(js*R k s)(* k a )) 


Now notice that 


$ s fc • M = 0 [J(V(j s * p s fc ))] 3 = 0 
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Hence 


• J($* - V(j s * R k a ))I < ||V 2 i?s | < |$J| max Ixl 

X£il2 


It follows that 


hence 


So we obtained 


d\®s(t,Vo)\ 

dt 


< max |x| t G [ks, (k + l)s] 


|*S((* + l)s,j/o)| < D k s + s • max \x\ < D k+1 


supp v k+l c B 2 D k +1 x [-£,-<5] 

Also we can define P k+1 , R k+1 , h k s +1 , 7 J+ 1 , w(? +1 , u (? +1 in the same way as above. 
Denote P s {t, -),R s (t, ■ ),h s (t , •), 7 s (t, •) to be above defined piecewise constant func¬ 
tion in time, i.e, P s (t, •) = P k t £ [fcs, (k + l)s), similar for others. 

Recall the definition of D, we have shown that if s < | 

supp v s (t,-) C B d+s ( 0) x 

Similar to the proof in [Bj , one has 

Wi(v a (t i,-).^(*2,*)) < (D + max|x|)||i/ 0 || i i( R 3)|ti -1 2 \ 

xen 2 

where is the 1-Wasserstein distance. Follow the proof of [ 8 j Theorem 5.3, 

we conclude that up to a subsequence 

v Sj —> v weakly in L r ([0,T] x R 3 ) 
v Sj {t, •) —>• u(t, ■) weakly in L r (A) 

6 L q ([0,T] x A) |KV)]|l«(A) = IN||l«(A) 


and 


Also 


Wi(y(ti,-),v(t 2 ,-)) < (£> + max |x|)||r' 0 ||L 1 (R3)|ti -t 2 \ 

xeii 2 

Since —> z/(t, •) narrowly as measures, we conclude by Theorem 2.25 that 

X7R Sj (t , •) — > VR(t, •) in L r ( A) and hence ->win L r ( A). 

UsjVsj (t, •) —> wz/(i, •) weakly in Z7(A;R 3 ) 

so we conclude and v satisfies the equation 

d t v + V • (wi/) = 0 ; i/(0, •) = vq 

in the sense of distribution, where w = J(y — Vi?(t, y)) and (P(t, ■),R{t , •)) is the 
minimizer of J^ t ^(P,R). The property of P,h follows from the stability result 
proved above and the narrow continuity of u t □ 
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2.2. Generalization to v with unbounded support. In this subsection, we 
will consider the case when v may have unbounded support and generalize the 
properties obtained in previous subsection. The result of this section will be used 
only in section 4, when the initial data has only L 2 , instead of L°° gradient. The 
ideas are quite similar, but certain complications arise. 

We will always take in this subsection A = R 2 x [— i, —5], and we assume v £ 
P 2 (R 3 ) with supp v C A. In this setting, E„(h, 7 ) is defined the same way as (2.1), 
J V (P,R) is defined as in (2.2). Obvious examples as P(x) = |( x 2 + x\),R(y) = 
\{yl + vl )shows 

(2.33) sup J V (P,R)> 0 

( p , r ) eAA 

Recall M v is the set of pairs (P, R) such that P(x) + R(y) > x ■ y with suitable 
integrability condition. See subsection 2.1. 

Suppose (P, R) £ M v , we can then use double convexification to define 


and 

we have 
Therefore 
then set Pq(x) 


R 0 {y) = sup (x-y- P(x)) 

X(zQoo 

Po(x) = sup (a; • y - Ro(y)) 

ye A 

Po < P and Rq < R 

J„(P,R) < J„(P 0 ,Ro) 
max(Po(a;), 2(ai 2 + x 2 )) and set 

Ri(y) = sup {x-y — Pq{x)) 


Pi{x) = sup(a: • y — Ri(y)) 
ye A 

Since P 0 (a:) > i(xf+x 2 ), it‘s easy to see R±{y) < \(y\ +y|). Since P 0 (x) + Ri(y) > 
x-y, we see Pi (a;) < P 0 (x) and the dehnition of Pi shows Pi (a;) + R\{y) > x ■ y. 

The proof of Lemma 5.1(H) in the appendix shows Pi (a;i, a: 2 ,0) > \(x\ +x 2 ). 
Since Ri(y) < R(y), we have 

j [\{yl + vi) - R(y)\dv(y) < J \^(y\ + vl) - Ri(y)\dv(y) 

Put ho(x) = hp 0 (x ), which is well defined since — | < d X3 P < —S, then we have 

inf [ [ij( x\+xl)-Po(x)}dx = [ [hxj+xl)-P 0 (x)]dx = [ [hxl+xl)-P 0 (x)]dx 
2 J Qho 2 J nhQ 2 

= inf / \k{x\ + * 2 ) - Po(x)]dx < inf f \^(x\ + x\) - Pi(a;)]dx 
h >oJn h 2 h ^°Jn h 2 

Since Pi (a:) < Pq(x), which gives J„(Po,Po) = Jv(Po,Ro) < J v (P\,Ri). Summa¬ 
rizing above discussion, we get the following lemma. 


Lemma 2.34. Let ( P,R) £ M v , then there exists (Pi, Pi) £ M v such that 

(*) MP,R)< Ju(Pi,Ri) 

(ii) Pi(xi,x 2 ,0) > \(x\ + xl) and R\(y) < \(yl + yl) 

(Hi) Pi ( x ), Pi (y) are convex conjugate over floo,A 
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Let ( P n ,R n ) be a maximizing sequence of J V (P,R). By the above lemma, one 
can assume ( P n ,R n ) has additional properties (■ ) above. Besides, we can 
also assume J v (P n , R n ) > 0 by (2.33). Next we will derive some estimates for the 
maximizing sequence, which allow us to pass to limit and prove the existence of at 
least one maximizer. 

Lemma 2.35. Let (. P n ,R n ) £ M v be a maximizing sequence of J U (P, R) with the 
properties (ii), (in) in previous lemma, Suppose J v (P n ,R n ) > 0 ,and put h n = hp n , 
then there exists a constant C = C^M^zv), c), f2 2 ), such that 

(2.36) \\hn\\ L 2 {Q2) ,\\Pn(;0)\\ L ^U 2 )<C 

and constant Ck = C(K , M 2 (zz), < 5 , fi 2 ), such that 

(2-37) I |-7n 11 l 2 (n 2 x [o,ic]) < C K 

We have also 

(2-38) ~max\x\\y\ - - Rn ^ - 3^1 + 

where Co has the same dependence as C and in particular H-RralUq*/) < C 
Proof. We start with the following estimate 

0< J v (P n ,R n ) = [\^(yl + yl)-R n (y)\dv(y)+ini /[i(x? +x 2 2 ) - P n (x)\a h (x)dx 
By taking h = jily, and noticing that — j < d X3 P n < — S , we see 


above < J\^(y{ + y\) - R n (y)\dv(y) + 

< ({^(yl + yi)-Rn{y)]du( y )+ f 

J ^ J fioxto. ... L . 


[-(x\ + x\) - P n (x)\dx 

1 Z 


l~(xl Px\) — P n (Xi,X 2 , 0) + -x 3 \dx 
2_i Z 0 


< ^ [M 2 (v)+ 2 max |af+ ^ £2 \n 2 ) }~ j R n(y)diz(y)—^ 


P n (xi,X2,0)dxidx2 


Now we notice 


P n (x 1,X 2 ,0) + Rn(y ) > Xiyi + x 2 V 2 > + A + vl + vl) 

Therefore integrating above inequality against xn 2 x dv, we obtain 
[ P n (xi, x 2 ,0)dx 1 dx 2 +C 2 (tt 2 ) [ R n (y)du(y) > -^M 2 (^)£ 2 (f2 2 )-^ 2 (f2 2 ) max |a;| 2 

Jfl 2 j * Z il 2 

So we deduce from (2.39) 

0 < i[M 2 M+2ma x Ixp+^T^l+ifMjfO+max P n (x u x 2 . 0) 

That is 
(2.40) 

[ P n (xi,x 2 , 0) < C (° 2 ^ [Af 2 (i/)+2max|a;l 2 + ^ J + ^ ^^ (M 2 (t/)+max |ar| 2 ) := C 0 

Jn 2 2 n 2 dL~(il 2 ) 2 n 2 
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Next we want to derive a bound for R n . By the L 1 bound derived above, we can 
find (a;", £ 2 ) € fi 2 , such that P n (xi, x %, 0 ) < c 2< $ 2 ) ,so 


^n(y) = sup (x-y- Pn(x)) > X™yi + X%y 2 - P„(x" .X 2 , 0) 


. 1 111 2C ° 

111191 -mthj 

To summarize 

2C 1 

~ m af |:C||2/| “ Z%|) - Rn{y) ~ 2 {vl+vl) 

So we have 

j\\{yl + y\) - Rn(y)\dv(y) < J\^{yl + yl) + max|x||y| + -^^y ]du(y) 

< M 2 (u) + max \x \ 2 + ° := Ci 

SA iP(Sl 2 ) 

Next we proceed to derive an L 2 bound for P„(-,0). Put h n = hp n . This is well 
defined since d X 3 P n < —5. Thus 

—C\ < inf f[\(xl + x\) - P n (x)\a h (x)dx = f [^-{x\ + x%) - P n {x)\<j hn {x)dx 
h>oj 2 J 2 

To proceed further, we notice that since P n (xi,x 2 ,h n (xi,x 2 )) = \{x\ + x 2 ), and 
for 0 < X 3 < h n . one has 

^{xj + x 2 ) < P n (x) - S(h n (xi,x 2 ) - x 3 ) 

So we conclude 

-Ci < -5 J [h n (xi,x 2 ) - x 3 ]a hn (x)dx = J h 2 n 


That is 

/ 

J fi 2 S 

Finally we notice 

2 ( x i + x 2 ) < P„( Xi, x 2 , 0) < — h n + P n (Xi, x 2 , h n ) = — h n + — (xf + x|) 

So we have ||P„(-, 0 )|| L 2 ( q 2) < C 3 . 

Since - \x 3 +P n {xi, x 2 , 0) < P n {x) < P n (xi,x 2 ,0)-6x 3 , we have \\P n \\L^(Q 2 x[o,K]) 
C{K , M 2 (u), fi 2 ) V n > 1 K > 0 □ 


The bounds derived in previous lennna allow us to pass to limit and prove the 
existence of at least one maximizer. By Lennna 5.5 in the appendix, we know 
P„(-, 0) is uniformly bounded on each compact set F of f l 2 and notice — | < d X 3 P n < 
—(5,P„are uniformly bounded on each set of the form F x [0,A']. Also we know 
they are equicontinuous on such a set and we can use Arzela-Ascoli to extract a 
subsequence which converges uniformly on each compact set of the form F x [0, A"]. 
Therefore, we can take a subsequence (not relabeled)P„, R n , h n , such that 
P n —» P uniformly on each compact subset of and in L r (Cl 2 x [0,AT]), for 
any r € [1, 2) and K > 0. 

h n —> h locally uniformly in fl 2 and in L r (Q 2 ), for any r € [1, 2) 

R n —> R locally uniformly in A and in L 1 (dv )(by dominated convergence) 
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Since for each n, we have P n {xi,X2, h n ) = ^(xf + x|), it‘s easy to see in the limit 
P(xi,X 2 ,h) = 2 ( x i + x|). This is by uniform convergence of h n ,P n on compact 
sets. So h = hp. 

We now need to show the limit (P, R ) is a maximizer. 

Lemma 2.41. Let P, R, h be as in previous paragraph, then we have 

J„(P, P)> lim sup J u {P n ,R n ) 

n—> oo 

in particular, (P, R) is a maximizer. Besides, 

P{x i,x 2 ,0) > ^(x? +^ 2 ) 

P(x) = sup(x • y - P(y)) 

2 /eA 

Proof. First since P n converges in L 1 {dv), we have 

J[\{yl + 2/1) - p«( 2 /)]^( 2 /) ->■ J[\{yl + 2/1) - P( 2 /)]<M 2 /) 

Next observe that on f2 An , we have P n {x) > ^(x^ + x|), so 0 < P n (x)ah n (x) —»• 
P(x)aj l (x) pointwise. So we can apply Fatou to see 


/ P{x)aj l {x) < lim inf / P n {x)<7h n (x)dx 

n—too J 

Combining the P 1 convergence of h n , we have 

/[^(x? + x\) — P{x)\aj l (x)dx > lim sup (\^-{x\ + x\) - P n (x)\a hn (x)dx 

J Z n—>• 00 J Z 

Therefore we see 

J U (P,R)> lim sup J v (P n ,R n ) 

n—>oo 

Since for each n, P n (x i,X 2 , 0 ) > \(x\ +X 2 ), one easily sees this is preserved in the 
limit. 

To see P(x) = sup ygA (x-j/ — R{y)), we first observe for each n,P n (x) = sup ygA (x- 
y — R n {y )) So P„(x) + Rn(y) > x • y, and the equality is achieved if y £ dP n (x). In 
the limit we have P(x) + R(y) > x • y. Therefore 

P(x) > sup(x • y - P(j/)) 

ye A 

To see the reverse inequality, fix x° € floo- Find a compact set K such that 
x° £ K C fioo. Now dP n (K ) must be uniformly bounded in A, by Lemma 5.5 in 
the appendix, so we can assume for some F compact dP n {K) C F C A, for any 
n. Given e > 0, by uniform convergence in compact sets we can take n = no, such 
that 

|P(x) - P„(x)| < e V x e K and \R n {y) - R{y)\ < e V y £ F 
Let p £ d P„(x°) C P, we then have 

P(x°) < P„(x°) + e = x° • p - R n (p) + e < x° ■ p - R(p) + 2e 


< sup(x° • y — R(y)) + 2e 
2/eA 


□ 
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On the other hand, if we define R{y) = sup^gQ (x ■ y — P{x)), then R{y) < 

R{y) < 3 ( 2/1 + 2 / 2)5 and so J U {P,R) < J„(P, P) and ( P,R ) is still a maximizer. 

Also because of the bound (2.36), (2.37) this bound will also be satisfied by P, so 
we can conclude R satisfies the bound (2.38). To summarize,we get 

Corollary 2.42. There exists a maximizer (Po,Po) of J V {P,R),such that 

(i) (P 0 ,R 0 ) convex conjugate over floo and A 
(m) P 0 (x 1 ,x 2 ,0) > \{x\ +x%) 

{in) Put h 0 = hp 0 ,then the bound in (2.36),(2.37),(2.38) hold true for P 0 , h 0 , Rq 
But the most useful properties will be the following 

(2.43) P{x) = sup (a: • y - R{y)) x £ fioo 

ye A 

and 

(2.44) P{xi,x 2 ,0)>^{xl+xl) 

The first condition ensures — i < d X3 P < —6 so that hp is a well-defined function 
on fl 2 . The second condition ensures P{x\,x 2 ,h) = ^{xf + x%). Besides they can 
be preserved in the limit. 

As before, we can prove the following property of maximizers. The argument is 
the same as Lemma 2.11, so we will omit the proof. 

Lemma 2.45. Let (P, R) £ AT be a maximizer of J„{P, R) and such that (2.43),(2.44) 
are satisfied. Let h{x) = hp{x), then VPjjcr/j = u. In particular, h = 1 

Now we consider E„{h, 7 ). Similar to what we proved in previous subsection, we 
have the following. 

Lemma 2.46. Let v £ P 2 (R 3 ) with supp v C R 2 x [— |, —5], then 
{i) E v {h, 7 ) > J U (P,R), for any {h, 7 ) G M and (P, P) G A f 

(ii) Suppose (Pq, Ro) be a maximizer of J V {P,R), and satisfies properties (2.43), (2.44), 
then ( i ) takes equality if and only if h = hp 0 and 7 = {id x VPo)j<T/, 

{Hi) E u {h,"f) has a unique minimizer (ho, 70 ) and there is a universal bound on 
h 

(2-47) ||ho||z, 2 (Q 2 ) < C(M 2 {v).6,Ll 2 ) 

{iv) If {Pq, Ro),{Pi, Ri) are two maximizers of J„{P, R), both satisfy (2.43),(2.44), 
then we have Po = Pl on f lh 0 = 0} 

Proof. The proof for Theorem 2.18 and Corollary 2.19 works here. The L 2 bound on 
h is a consequence of Lemma 2.35. Also recall that by Corollary 2.36 (ii)whenever 
h{xi,x 2 ) = 0 for some ( 27 , x 2 ) G fi 2 , we have exactly P(cci, a; 2 , 0) = \{x\ +x 2 )- □ 

Next we study the stability property of these optimizers under narrow conver¬ 
gence of v which remain bounded in P 2 (Ra). 

Suppose u”, v G P 2 (R 3 ), supp v n , v C A, and v n —» v narrowly with sup n>1 M 2 {v n ) < 
00 . Let {h n , 7 ”) be the minimizer of E„n{h, 7 ), {P n ,R n ) a maximizer of J„n{P,R) 
given by Corollary 2.42. Then sup n>1 W 2 {v n ,v) < 00 . Now noticing {Hi) of Corol¬ 
lary 2.36, , we can take a subsequence such that 

P n —> P uniformly on each compact subset of floo and in L r {il 2 x [0, A]), for 
any r G [1, 2) and K > 0 

h n —> h locally uniformly in f l 2 and in i r (fl 2 ),for any r G [1, 2) 
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P" —> R locally uniformly in A 

Put 7 ” = (id x VP n )$ahn and 7 = (id x VP)jjcr^, then above convergence implies 
7 " —>• 7 narrowly. Indeed for any g G Cb (R 3 x R 3 ), we have 

J 9 (x,y)d'y n (x,y) = j g(x, VP n (x))a h n(x)dx ->• J g(x,VP(x))a- h (x)dx 


= J g(x,y)d ; y 

The above convergence is due to L 1 convergence of h n and point-wise convergence 
of VP” 

Similar as before, the following semi-continuity result holds. 

Lemma 2.48. Let h n , 7 n ,/i, 7 6 e defined as in previous paragraph, then we have 

E„(h,j) < lim inf P y n(/i”, 7 ”) 

n—> 00 

Proof. This follows from the narrow convergence of 7 " since the integrand is non¬ 
negative. □ 

In the following lemma, we show that the limit (h, 7 ) and (P, R) obtained above 
are indeed optimizers of E v (h, 7 ) and J U (P,R) respectively. 


Lemma 2.49. The following statements are true. 

(i) P(x) = swp yeA (x ■ y - R(y)), and P(xi,x 2 ,0) > \(x\ + x\) 

(ii) (ft., 7 ) is the unique minimizer of E u (h,"/), 

(in) (P, R) is a maximizer of J„ (P, P). 


Proof. We first observe P(x) = sup yeA (x • y — R(y))- The argument is the same as 
Lemma 2.41. 

To see (h, 7 ) is the minimizer, and (P, R) is a maximizer, we prove it by showing 
E v (h, 7 ) = J„(P, P). Indeed, for each n, 

(2.50) P"(x)+P”(VP”(x)) =x-VP"(x) ^a.exefloo 
Passing both sides of (2.50) to limit, we then have 

(2.51) P(x) + P(VP(x)) = x ■ VP(x) C 3 a.e x G floo 
This means 

(2.52) P(x) + R(y) =x-y 7 a.e (x, y) G x A 
Therefore 

K(h, 7 ) = J [^(xl + xl + yl + yl) - x ■ y\d*f(x, y) 

= J [§( 2 /i + 2/2) - #( 2 /)]<M 2 /) + J [^(*1 + ^2) - P(®)]^ (x)dx 
Since h = hp 

J 1 + X2) - P(z)K(x)dx = inf J \^(x\ + x^) - P(x)]tr/ l (x)dx 


We obtain E„(h, 7 ) = J„(P, P) 


□ 


We are now ready to prove the following stability result. 
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Theorem 2.53. Let v n ,v G P 2 (R 3 ), suppv n ,u c A, with v n —> v narrowly 
and sup n>1 M 2 (v n ) < oo. Let (h n , 'y n ),(h, 7 ) be the unique minimizer of E u (h,"f), 
E v (h,j) respectively, let (P n , R n ),(P , R) be a maximizer of J v n(P,R), J V (P,R) 
respectively which satisfies (2.43),(2.44), then the following convergence are true: 

(i) h n —>• h in L r (Q 2 ), for any r G [1, 2). 

(ii) 7 " —»• 7 narrowly. 

(Hi) £(P", VP n )a h n -a £(P, VP)^ in L 1 ^), for any £ G C b (R x R 3 ) 

If one further assmes W 2 (u n , v) —> 0, then (i) can be improved to 
(iv)k n h G L 2 (fl 2 ) 

Proof. We can see from previous lemma that for any subsequence of (h n , 7 ™), there 
is a further subsequence, say (h nj , j Hj ), such that h nj converges in L r (Ll 2 ), for 
any r G [1,2), 7 "^ converges narrowly, and the limit is the unique minimizer of 
E„(h, 7 ). This is sufficient to conclude the whole sequence must converge. This 
proves (i), (ii). 

The argument for (Hi) is similar. We first show any sequence has a further 
subsequence which converges C 3 — a.e to £(P, VP)^. Indeed, let ( P n ,R n ) be 
maximizers of J„n(P, P) given by Corollary 2.42. Let £(P", VP")<T/jn be a subse- 
quence(not relabeled), by what has been discussed, we can take a further subse¬ 
quence, say (P n '>, R n i) which converges locally uniformly on ,A respectively to a 
maximizer (P,P), this maximizer will satisfy (2.43),(2.44), also £(P" 3 , VP" 3 )<j/ ( >q 
will converge C 3 - a.e to f(P,X7P)a ~ h , but £(P",VP"H» = £(P", VP> /l », 
£(P, S7P)aj l = £(P,VP)a~ h because of the uniqueness property proved in lennna 
2.46 (iv). 

To see such a convergence is in P 1 , we just need to observe that sup n>1 f^ X3>K y <?h n dx 
0 as K —} 00 . Indeed 

[ cr h ndx = f (h n - I<) + dx 1 dx 2 < - ( ° 2) 

J{x 3 >k} Jq 2 k 

Recall the universal bound on ||/i n Hz, 2 (n 2 ) asserted in (2.47), we proved (Hi). 

Now we assume W 2 (v n , v) —>• 0. Let (P n ,R n ) be the above chosen maximizers 
of Ji, n (P, R) given by Corollary 2.42 and nj is the above chosen subsequence which 
converges uniformly on compact sets. Because of the bound (2.38), and the assumed 
convergence W 2 (u n , v) —> 0, we have 

J\\(yl + 2 / 2 ) - (v)\dv n (y) ->■ J [^( y\ + vl) - R(y)\dv(y) 

while by Fatou(notice the integrand is nonpositive) 

/ (m? + * 2 ) — P( x )\ a h (x)dx > lim sup [ [—(ar? + x 2 ) — P n (x)\crh n (x)dx 

J 2 n—too J 2 

Hence 

J V (P,R) > lim sup J v n(P n ,R n ) 

n—> 00 

Recall Lemma 2.46(f), (ii) and Lennna 2.48, we actually have under the stronger 
convergence of u n 

E v (h, 7 ) = lim E„n(h n , 7 ") 

n—> 00 
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As a result of this,and noticing that \x\yi + x^y^i | < + 2 /f),we have 

J(-y 3 )x 3 dp n {x,y) -> J (- V 3 )x 3 (J~ h d ; y{x,y ) 

Recall 0 < Ar 3 < (— y 3 )xo,we obtain 

J x 3 <Jh n {x)dx —»• J xz<r- h (x)dx 

which implies 


n\2 


W 2 


□ 


/ (*") 

</ f 22 J 

Combined with (i), we get L 2 convergence of h n 

To conclude this section, we observe that for the maximizers found for J(p(P, R) 
in subsection 2.1 are also maximizers for J„(P, P), upon suitable extension. This 
fact is contained in the following lemma. 

Lemma 2.54. Let v £ P(R 3 ) with supp vcA, where A is as in (2.26). Let H 
he chosen as (2.32). Let (P, R) he the unique pair of maximizer given by Corollary 
2.20. Define P(x) = sup ygA (a; • y — R(y)) x £ to be the extension of P to fioo. 
Then 

(i) (P, R) are convex conjugate over A 

(ii) (P, R) is a maximizer for J„(P,R) 

Proof. To prove (*), we just need to check 

R(y) = sup (x-y- P(x)) 

This is routine. To see (P, R) is a maximizer for J„(P, P), we recall since (P, R) 
is assumed to be a maximizer of j(f(P,R), inf E v (y, h) = j(f(P,R) by Theorem 
2.18. But we have j(f(P,R) = J„(P,P) because hp = hp < H as guaranteed by 
the choice of sufficiently large H. □ 

3. Existence of Lagrangian solutions 

In this section, we prove the existence of weak Lagrangian solutions when the 
initial dual density is absolutely continuous, using the properties of geostrphic func¬ 
tional E v {h, n f) already proved in subsection 2.1. The proof here is similar to [5]. 

3.1. Basic definitions and main result. First let's define the notion of admis¬ 
sible initial data and weak Lagrangian solution. 

Fix H Ql 5 > 0, 

Definition 3.1. Given ( ho,Po ) with ho £ W 1,0O (fi 2 ) and Po £ W 1, °°(LLh 0 ) with 
Hq > maXxgrj 2 Po(x 1 .x 2 ,o) , sa y it j s an admissible initial data if the following holds 
true. 

(«)0 < h 0 < H 0 , f n h 0 dxidx2 = 1, 

(ii) Po is convex, and dPo(Ll 2 x [0, Po]) C R 2 x [—|, —<5] and bounded. 

(iii) P(xi,X 2 ,h 0 (xi,X 2 )) = \{x\ +x\). 

Remark 3.2. Above assumption H 0 > maX:rgsi 2 Po(x 1 .x 2 , 0 ) ac t; Ua ny guarantees h 0 < 

H 0 
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The above definition guarantees at least on f lh 0 , Po is the restriction of some 
maximizer of J Vo (P,R), where no = VPoftOAo. This is shown by the following 
proposition. 

Proposition 3.3. Let Do = ||VPo||loo + 1. let Ao = Bd 0 x [— j-,— 5] , H as in 
(2.32). Let (h 0 ,P 0 ) be admissible initial data. Let vq = VPojcr; lo so that supp vq C 
Ao, Let (Po,Po) be the unique maximizer of J^(P, R) which satisfies the properties 
(i), (ii) in Corollary 2.20, (ho, 7 ) be the unique minimizer of E vo ( h , 7 ), then Pq = Po 
on Cl ho U{x 3 = 0}, and h 0 = h 0 ,7 = (id x VPo)(,a ho 


Proof. Define 

Then one has 


R 0 (y) = sup (x ■ y - P 0 (x)) 
x€Qh 

P 0 (x) + R 0 (y) >x-y x e il H y G A 


Also observe that 

h 0 (xi,x 2 ) = hp 0 (x i,x 2 ) 

By (Hi) of above definition. 

We just need to show E Vo (ho,Jo) = J„ 0 (Po,Ro) 
Indeed, by definition 


J„ 0 (Po,Ro) — 


h(yi+y2)-Ro(y)\My)dy+ inf 

Z 0 <h<H 


r l 


(x\+xl)-Po(x)\a h (x)dx 


Jn 


= [ \k(vl + 2/1) - Po(y)\My)dy + f \\(x\ + x%) - P 0 (x)\a ho (x)dx 
JA 1 Jn^ 4 

= [ i\(xi+xl)-Po(x) + \(yl + yl)-Ro(y)\dio(x,y) 

JflooXA 1 1 

The second equality is due to above observation, and the third equality is because 
70 has <Jh 0 and uo as marginals. 

We will have shown j(f o (P 0 ,Ro) = E Vo (ho,Jo) provided we can show P 0 (x) + 
Ro(y) = x ■ y 70 a.e 

Indeed, let's show for xo 6 ilh 0 , such that Po is differentiable at Xo, we have 

P 0 (x 0 ) + P 0 (VP 0 (x 0 )) = x 0 • VP 0 (x 0 ) 


This is implied by 


Po(x) > Po(xo) + VP 0 (x 0 )(x - X 0 ) X e floo 

Since Po is convex in floo. So we have shown j)f(Po,Ro) = E VQ (ho,Jo)- So (ho, Jo) 
defined above is the unique minimizer. 

Now let (P 0 ,Po) be the unique maximizer of jff 0 (P,R ) given by corollary 2 . 20 , 
then we have 

7o = (id x VPo^Uh,, = (id x 37P 0 )ta ho 

and 

ho = h^ 0 = h\ 

The same argument as in Corollary 2.19 implies Po = Po on fb , 0 . Also for (xi, x 2 ) G 
il 2 such that ho(x\,x 2 ) = 0, we have P(x’i,X 2 , 0 ) = \(x\ +x|) 

□ 


Next we define Eulerian solutions for the system in physical space. 
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Definition 3.4. Let T > 0 ,H > 0, Let P(t,x ) £ L°°([0,T],W 1,oc, (^h)), such that 

1 1 —> P(t, •) is convex Vt € [0, T] 

Let u £ L 1 ([0, T] x Cl H ). 

Let ft € L°°([0,T], W 1,00 (fi 2 ))be such that 


0 < ft < H / ft = 1 


Then we say the triple (P, u, ft) is a weak Eulerian solution if the following is 
satisfied. 

(i) 


/ / u • V</>(&, t)dxdt + / / 

*^0 *^0 

V<j> € ct([0,T) X n H ) 


d t (j)(t, x)dxdt + / 0 (O,a:)da; = O 


(■ ii)P(t , &i, x 2 , ft) = 5(^1 + 

(Hi) 

n S7P(t, x)d t ip(t,x)dxdt + f 
!h(t-to 

/•T 


VP(f, a;) • (u • V)ip(t, x)dxdt 


[ [ [J(yP(t,x) — x)] ■ ip(t,x)dxdt + f VPo(a;)^(0, x)dx = 0 

J 0 J£lu(+ .'i 


/O 

WeCl({0,T)xn H :R 3 ). 


In the following we show that a weak Eulerian solution with sufficient regularity 
gives a classical solution. This justifies our definition. 

Proposition 3.5. Supposeu £ C 1 ([0, T]xf2#), ft £ C 1 ([0,T , ]xf! 2 ), P £ C ,2 ([0,T]x 
fl#) and (P,u,h)is a weak Eulerian solution, then they solve the equation in the 
classical sense. 


Proof. First we wish to deduce from (i) above the divergence free of u and the free 
boundary condition. Indeed under our regularity assumption 

/ / u • x)dxdtdxdt 

J 0 J Q h 


'0 j dQ.h-{x 3 =h} 


(u • n )<j>(t, x)dSdt + 


>0 J {x 3 =h} 


(u • n )<f>(t, x)dSdt 


Jo Jci h 

The second term above can be written as 


(V • u )4>(t, x)dxdt 


/ / (u-n )4>(t,x)dSdt = / / (~uid Xl h - u 2 d X2 h + u 3 )(j)(x,t)\ X3 =hdxdt 

JO Jfx*=h\ Jo Jflo. 


'0 J {xz=h} 

For the rest of the terms 


II dt4>(t, x)dxdtdx 
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n l>h(t,X 1 , 2 : 2 ) pT n 

(d t / <j)dx 3 )dt - / d t h<l)\ X3 =hdx 1 dx2dt 

i 2 Jo Jo Jn 2 


Combining terms to obtain 

pT r 


(\)dx —I / dth(j)\ X3 =hdxidx 2 dt 


/o Jn 2 


/ 0 J dClh — {x3=h} 

pT 


(u • n)(/>(t, x)dSdt — / / (V • u)0(t, x)dxdt 


/0 J £ih 


+ / (~d t h - uid xi - u 2 d X2 h + u 3 )<j>(x, t)\ X3=h dxdt = 0 

0 J Q 2 

By choosing appriorate test funtions we see that 

V • u = 0 in Vlh 

As well as 

u • n = 0 on dn h - {x 3 = h} 

and 

dth + u\d Xl h + u 2 d X2 h = u 3 on {.T 3 = h} 

Finally we recover the equation D t (VP) = J(VP — x) 

First observe that 


/ / VP(t,x)d t ip(t,x)dxdt 

Jo JQh 

T ph(t,x 1,2:2) ph(t,x 1,^2) 

VP(t,x)ip(t,x)dx 3— / d t (VP)(t,x)ip(t,x)dx 3) 

Jo 

T 


= / dx\dx 2 / («9 f / 

Jn 2 ./o Jo 

- / dx\dx 2 / \7P(t,x)il)(t,x)d t h(t,x)\ X3 =h 
J Qo J 0 




VPqOzOVKOj — / / dt(VP)(t,x)'i/)(t,x)dxdt— 


/o 


dxidx 2 / VP(t, x)d t h(t , x)| X3= h 
Jo 


while the term 


VF(t,i)(u • V)ip(t,x)dxdt 


>0 JOh 


= f f V ■ [(VP ■ i/))u\ — (u-V)(VP(t,x)) ■ ip(t,x)dxdt 
Jo Jn h 

= I f (VP ■ ip)u ■ ndSdt — f f (u • V)(VP(t, x)) ■ ip(t, x)dxdt 

Jo J d£lh, Jo J £lh, 

= f [ (-uid Xl h-u 2 d X2 h+u 3 )(VP-ip)(t,x)\ X3 = h dxdt- f ( (u ■V)(VP(t,x))-ip(t,x)dxdt 
Jo Jn 2 Jo Jn h 

In the above ,we used divergence free of u as well as the boundary condition u • n = 

0 on dflh — {^3 = h} 

Now collect terms and use the free boundary condition to see the resulting equa¬ 
tion is exactly what we want. □ 
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Now we define the notion of weak Lagrangian solution with admissible initial 
data. 

Definition 3.6. Let T > 0 ,q £ (1, oo), Let h be such that 

h(t, Xl ,x 2 ) £ L oo ([0,T),W 1 ' oo (n 2 ))f|C°([0,T] x n 2 ) 

Let P(t,xi,x 2 ,X 3 ) defined on Qh be such that 

1 1 —> P{t , •) is convex Vf £ [0, T\ 

and additionally 

P £ L°°([0, T], W 1 ’°°(fi H ) f| C([0, T] x W x ’ r (n H )) 

which assumes initial data (ho,Po) on f \ 0 1 J {®3 = 0 }, i.e h(0,x) = h 0 (x) x £ fl 2 
and P(0,x) = P 0 ( x) x £ fl ho |J{;r 3 = 0} 

Let F : [0, T] x f \ 0 —> fin be a Borel map and such that for some 1 < r < oo. 

F £C([0,T]-L r (o ho dC 3 , R 3 )) 

Then we say the triple (h, P 1 F) is a weak Lagrangian solution with initial data(/io, Po) (admissible 
in the sense above) if the following holds: 

(i)For each t £ [0,T], 



+ / Z(0, x)<j)(0, x)dx = 0 

J Qh 0 

V<j>(t,x)£Cl([0,T)xn H ) 

Remark 3.7. (f)We only assume P achieves the initial data on f2/j 0 |J{a :3 = 0} 
instead of fl# because only the region £lh 0 is physically relevant. 

(n)In the ( v ) of above definition, it is possible to choose test function cf> such 
that d t (j) £ L°°([0,T) x f Ih) without assuming V x 0 exists. Indeed, we may define 
(f>(—t,x) = ^(0, x) for t > 0 and convolve with J s {t,x) with s small such that 
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(j) s : = (j)*j s is a legitimate test function. Then dt<t> s ,<t>s converges a.e in [0,T] x floo 
to dt(j>, (f> as s 0 then the result follows from dominated convergence. 

(ra)In the ( iv ) above, note that 

U {<} x Qh(t,-) = G [0, T] x fl H \0 <x 3 < h(t,xi,x 2 )} 

te[o,T] 

Of course we must show above definition does not lose any information, i.e we 
need to show with additional regularity assumption, weak Lagrangian solution gives 
weak Eulerian solution. 


Proposition 3.8. Suppose (h, P, F) is a weak Lagrangian solution, suppose also 
that dtFt G L°°([0,T] x Qh), Define u (t,x) = dtF t (Ff(x)), then (P, u, h) gives a 
weak Eulerian solution. 


Proof. We only need to check that (i) and (Hi) of the definition of Eulerian solution 
are satisfied. 

First notice that F t ^ah 0 = &h- Thus for <f G C).([0,T) x Ll^), we have 


f ( u • V<j)(x,t)dxdt = f f d t F t (Ff(x))'\74>(x,t)dxdt 
Jo Jn h Jo J 

d t F t (x)'\7(j)(Ft(x), t)dxdt 
d t (J>(F t (x),t)) - d t (j>(F t (x),t)dxdt 


n 

>0 Jn hn 


f 

1 o Jn hn 


/ <j>(x, 0)dx — / / dt<f>(F t (x),t)dxdt 

'fi h n Jo J^h a 


(j>(x, 0 )dx — 


/o Jn h 


dtf>(x, t)dxdt. 


This verifies (i). 

In the second line above, we used F t ftOh 0 = <Xh(t,-) an d also Ft* (Ft (a;)) = x Oh 0 a.e 
In the third line above, we used 


for a ho a.e x, d t (fi(F t (x),t)) = d t <t>(F t (x),t) + V</>(F t (x), t)d t F t ( x) 
since t \—> F t (x) is Lipscliitz. 


Now we verify (Hi). Indeed 

[ [ VP(t, x)d t 4>(t, x)dxdt = [ 


Z(t, x)dt'fi(t, F t (x))dxdt 


Z(t,x)[d t (if(t,F t (x)) - V'lfft,F t (x))d t F t (x)]dxdt 


Z(t, x)dt(i/>(t, F t (x)))dxdt — 


(u • S7)ip(t, x)S7P(t, x)dxdt 


In the second line above, we used the fact that since d t F t (x) G P°°([0,r] x 
^H),d t (ip(t, F t (x))) G P°°([0,T] x fiff) and usual chain rule holds. 
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Now we choose test function as ip(t,Ft{x)) in the definition^) above to get the 
following. This is justified because of the remark after the definition above. 

f f Z(t,x)dt{il>(t,F t (x))dxdt 
Jo Jn hQ 


'o JQh 0 
r T 


J(Z(t,x) — F(t,x))ip(t, F t (x))dxdt — / Z(0, x)if(0, x)dx 

" ^/lg 

J(VP(t, a:) — x)ip(t, x)dxdt — / \7Po{x)ip{0,x)dx 

j Qhr, 


Jo Jn h 

Put things together, we get {Hi). 

Now we can state the existence result of weak Lagrangian solutions. 


□ 


Theorem 3.9. Let T > 0,,1 < q < oo, and admissible initial data ( ho,Po ) be 
given, suppose also 

is 0 ■■= VPotJCTfto e L q (R 3 ) 

Suppose also that 

(3 - 10) H > ZW + 2dia 7 ^ [l l V *H- r + m nf l*K T + 2) + 2] 

then there exists a weak Lagrangian solution {h, P , F) on [0, T] x CIh- Moreover the 
function Z{-,x ) G W 1,00 (R 3 ) /or a.e a: G f2^ 0 and the equations are satisfied in the 
following sense: 

d t Z{t , x) = J{Z{t , x) — F(t, x)) for £ 4 — a.e in {t, x) G [0, T] x Llh 0 

Z(0,x) = \7 Pq{x) for C 3 — a.e in x G flh 0 

3.2. Lagrangian flow in dual space. Next we study the Lagrangian flow in dual 
space. It is similar to [5] section 2.3 
Recall that is 0 = VPo^ah 0 ,so 


(3.11) 

supp vo C -Bd o ( 0) x -£] 

where 


(3.12) 

L>o = ||VP 0 ||loo + 1 

Put 


(3.13) 

D = Dq + T max |x| + 1 

and define 

A = B d (0) x 

(3.14) 


Choose FI as in (3.10). Let h, P , R, v be the dual space solution given by Theorem 
2.30 where we have chosen the parameters T, H, D, D 0 as here, recall that is{t, •) 
satisfies the transport equation 

d t v + w • S7v = 0 i/(0, •) = Vo 

where w = J(y — \7R{t,y)) is divergence free. 

Then it follows from Theorem 2.30 that 

w G L“ c ([0, T\ x R 3 ) w G L°°([0, T],BV{B(0,R))) VP > 0 
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Here we naturally extend w to R 3 by the same formula above. 

Since v is supported in Bd( 0) x [—y, — <5] x [0, T] we can modify w outside Rd(0) x 
[— y, —5] such that the modified w G L°°([ 0, T] xR 3 ) and w G L°°([ 0, T],BV(B( 0, R))) \/R > 
0 and V • (w (t, •)) = 0. To construct w, let £ £ CP^R 3 ) be a cut-off function such 
that 

C(s) = 1 for |s| < D + - £(s) = 0 for |s| > D + 1 + 1 0 < £ < 1 
o o 

Then define 

H{y) = (C(|j/i|)2/i> C(lz/2 D2/25 C(lz/ 3 D2/3) 

Then we take 

w = J(H(y) - X7R(t,y)) y G R 3 

Now we can apply Ambrosio‘s theory pQ to w to get a Lagrangian flow $ in dual 
space and establish the following lemma. 

Lemma 3.15. Let w be defined as above. Then there exists a unique locally bounded 
Borel measurable map $ : [0, T] x R 3 — > R 3 , satisfying 
(*)$(•, y) G W 1,oo ([0,T];R 3 ) for a.e ye R 3 
(m)$( 0 , y) = y for C?-a.e y G R 3 

( iii) For a.e (■ t,y ) G R 3 x (0, T) 

dt$(t,y) = w (t,$(t,y)) 

(iv) d>(t, ■) : R 3 —>■ R 3 is a C 3 measure preserving map for every t G [0,T] 

Remark 3.16. Notice that by our definition of w, we have w 3 = 0. By (i) and (Hi) 
above, one sees that 

* 3 (t,y) = y Vi G [0, T] 

for C 3 — a.e y 

Lemma 3.17. Let w be defined as above, and let $ be the flow in previous lemma, 
then 
(*) 

Hfy) C B 2 d ( 0)) x [-1-6] for a.e (; t,y ) G [0,T] x VP 0 (^oo) 

In particular, 

dt^(t, y) = w(t, <f>(i, y)) for a.e (t, y) G [0, T] x VP 0 (H oo ) 

(■ ii) There exists a Borel map : [0, T]xR 3 -> R 3 such that for every t G [0, T] 
the map : R 3 —> R 3 is C? measure preserving, and such that o$ t (y) = y and 
o <f> t * (y) = y for a.e y £ R 3 

The proof of this lemma is almost identical to the lemma 2.9 om , so we omit 
the proof. 

Proposition 3.18. Let Qn,T,q be as in Theorem 2.30. Let h,P,R,v be the dual 
space solution obtained in that theorem. Let w, $ be as in above definition. Let 
D = ||VPo||l~ + (T + l)maxfi 2 |x| + 1, then for t G [0,T] 

v(t, ■) = 4> t |juo 

Moreover, for t G [0, T], we have 

v(t,y) = 1 'o($t(y)) a.e y 
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The proof is again similar to Proposition 2.110. Keep in mind that we have 
strong convergence of the transporting vector w s . to w according to the proof of 
Theorem 2.30. This enables us to use the stability result Theorem 5.5 of [T] 


3.3. Lagrangian flow in physical space. We want to define a Lagrangian flow 
in the physical space F : [0,T] x Q/, 0 —>• fl# by defining F t : f l ho —> fl H by 

F t := VR t o$ ( o VP 0 

Of course, we need to check above formula is well defined. 

Lemma 3.19. The right hand side above is defined Cfi — a.e. For any t £ [0, T], 
above is defined £ 3 — a.e. and F : [0,T] x Qh 0 —> Qh is a Borel map. 

This is checked in the same way as Lemma 2.12 of [5]. So we omit the proof. 

In the following, we collect some results which can be proved in the same way as 
[5j. they correspond to Proposition 2.13-2.16 in that paper. 


Lemma 3.20. 


F{ 0, x) = x a.e x £ f \ 0 
Lemma 3.21. For every t £ [0, T], we have 

F n <Jh 0 = °7i(t,-) 

Proposition 3.22. For any to £ [0 ,T\,any r £ [l,oo), one has 


lim [ \F t {x) — F to {x)\ r dx = 0 

t-fto,t6[0,T] Jci hQ 

Proposition 3.23. There exists a Borel map 


F* '■ U W x ^h 

te[o,T] 

such that for any t £ [0, T] 

F t (F*(x)) = x <r h ( t> .) - a.e and Ff(F t (x)) = x a ho - a.e 


Proposition 3.24. Put Z(t,x) = 
Z(t, x)df (j)(t, x)dxdt 


fl 

' 0 J Qh, 


o 


' VPt(F t (x )), then Z(t, x) satisfies 

+ f ( J(Z(t,x) — F(t,x))(f>(t,x)dxdt 
Jo J^h 0 


+ / Z(0, x)(f>(0, x)dx = 0 

V0(t,®)€C , c 1 ([o,r)x n H ) 

Moreover, possibly after chanqinq Z(t,x) on a neqliqible subset of [0,T1 x Qj ln , we 
have ZM £ w‘ '“([0,r];R.’) for a.e-x e SlJ and 

d t Z{t, x) = J(Z(t, x) — F{t , x)) for £ 4 — a.e in ( t , x) £ [0, T] x flh 0 
Z(0,x) = VPq(x) for C? — a.e in x £ Qh 0 


Proof. By Lemma 3.20, we have 

Z(t, x) = VP t o F t (x) =$ f o VP 0 (i) 

Except for a £ 4 negligible set. 

So after redefining Z on a negligible set, we may redefine Z{t,x) = 4> t o VPq(x) 
and we will prove this version of Z has all the properties claimed. 
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By Lemma 3.16, we know that $(-,y) G VL 1 ,oo ([0, T]; R 3 ) for a.ey. Since 
S7Po(,<Jh 0 « £ 3 , one has 


t^$ f (VP 0 W) G W 1,oo ([0,r];R 3 ) for C 3 - a.e x G fi ho 

Let 1V 3 be such that C 3 (N 3 ) = 0 and for y G N 3 ,(i),( ii) of Lemma 3.16 holds, and 
such that for a.e t G [0, T ], y) = w(t, y)) holds. Then for such y, one has 

y) = y + [ w(s, $(s, y))ds Vi G [0, T] 

Jo 

Due to the same reason as above and note that |VPo(x)| < D Vx G £lh 0 by our 
choice of D, one can conclude for C 3 — a.e x G f lh 0 


$(f,VP 0 (x)) 


VP 0 (x) + f w(s,$(s,VP 0 (x))ds Vi G [0,T] 
Jo 


Therefore 

Z(0,x) = x C A — a.e x G fbi 0 

Also for such x, one lias 


d t$(f, VP 0 (x)) = w(t, $(t, VP 0 (x))) 


= VP 0 (x)) - VP t o o VP 0 (x)) 

= J(Z(t, x) — F(t, x)) a.e f G [0, T] 

So we have shown now 

d t Z{t, x) = J(Z(t , x) — F(f, x)) for £ 4 — a.e in (i, x) G [0, T] x 

Next we notice since we assumed Pq to be Lipschitz on £Ih- VPo(fiij) will remain 
in a bounded set. And the flow map $ is locally bounded on R 3 x [0, T], therefore 
Z(t,x) =$ t o VPo(x) will be bounded for fin x [0,T]. So we can multiply a test 
function to above a.e identity and integrate by parts in t to get the distributional 
identity. □ 


Proof of theorem 3.9: 

Let D 0 = ||VP 0 ||loo + 1, set A 0 = B Dq ( 0) x [- 55 ,-f], and put u 0 := VP ot jtT^ 0 . 
Then supp u 0 C A 0 . Also put D = D 0 + (T + 1) maxn 2 |x| + 1 and A = Br>{ 0) x 

Choose H as stated in Theorem 2.3. Then such H also satisfies 


H > 


2 | 2 maxo 2 |x|+ 2 P 

£ 2 {Sl 2 ) + 6 


• diam 


as required by Theorem 2.30. 

Now let (h, P, R) be the dual space solution corresponding to the initial data u 0 
given by Theorem 2.30, and let F be as defined as in the begining of this subsection. 
Then (h, P, F) gives a weak Lagrangian solution with initial data (ho, Po). 

Indeed, that (h,P) assumes initial data (ho,Po) on f lh 0 = 0} is guaranteed 

by Proposition 3.3. property (i),(ii) comes from our construction of dual space 
solution, other properties follows from the lemmas and propositions listed in this 
subsection. 
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4. Existence of relaxed Lagrangian solution with general initial 

data(VPq e L 2 ) 


4.1. Definition of generalized data and main result. In this section, we prove 
the existence of relaxed Lagrangian solution in a similar way as was done in m- 
The relaxed Lagrangian solution is an even weaker notion than weak Lagrangian 
solutions defined in previous sections, but it will allow for more general initial data 
(ho,Po), in particular, we will no longer require VPojct/jo << £ 3 or have compact 
support (still we assume supp C { —j < X 3 < —5}). To motivate the definition, 
recall that the weak Lagrangian solution given by Theorem 3.10 will satisfy the 
additional property 

dtZ(t, x ) = J(Z(t , x) — F(t, x)) for £ 4 — a.e in ( t, x) £ [0, T ] x flh 0 

Z(0,x) = \7 Pq(x) for £ 3 — a.e in x £ flh 0 
Recalling the definition of Z(t,x), this implies for £ £ C ,1 (R 3 ), we have 
dtt{VPt{F t (x)) = V&VP t (F t (x)yj(yP t {F t (x))-F(t,x)) for £ 4 -a.e in (t, x) £ [ 0 , T]xfl ho 


and has initial data £(VPo(x)). 

Thus if we define a Borel family of measures [0, T] 3 t 1 —> da t (x,x) := (id x 
Ft)$<Xh 0 and put da = datdt, then at least formally we can obtain from above 
almost everywhere defined equality that 
(4.1) 

r T r [>T /* 


/ / €(VP t (x))d t il>(t,x)dat(x,x)dt+ / Vf(\7P t (x))-J(\7P t (x)-x)ip(t,x)da t (x,x)dt 

J 0 •/ n't. J 0 J ns 


+ [ £(VPo(a;))V’(0, x)dx = 0 V£ G C^R 3 ) 

•'Qh.Q 


Above discussion motivates the following definition of relaxed Lagrangian solutions. 


Definition 4.2. (h 0 ,P 0 ) be admissible initial data in the sense of Definition 4.3 
below. Consider a Borel function P : fx [0,T] —> R, such that P(t, •) is convex 
for any t £ [0, T] and h : VL 2 x [0, T] —> R such that h(t , •) £ L 2 ( ^ 2 ) and a family 
of Borel measures [0,T] 3 t 1 —> a t (x,x) £ P(Doo x Rqo)- Let a be given by 
da = datdt. We say that the triple (P, h , a) is a relaxed Lagrangian solution if the 
following holds 

(i) h > 0 f n2 hdx\dxi = 1 V<€ [0,T] 

(ii) P t (xi,X 2 ,h) = \(x\ +x\)(S7Pt £ L 2 (n ht ), and < d X 3 P t < -5 for any 
t € [ 0 ,T] 

(iii) proj x a t = o ho proj x a t = cr h ( t ,-)> «o = (id x id)^a ho 

(w)(4.1) above holds true. 

(v)P(0, •) = P 0 on fl ho h( 0, •) = h 0 

Since we derived (4.1) only formally, we need to check (4.1) makes sense. 


Lemma 4.3. The left hand side of (4-1) is well-defined for any da(t,x,x) = 
dat(x,x)dt with at satisfying (Hi) in Definition 4-2 . 


Proof. The main issue comes from the fact that for each fixed t, VPj is not an 
honest function but is defined only a.e, We have to check that different choice 
of Borel representative of VP* does not affect the integral in the left hand side. 
Actually we will show for each fixed time slice t, the inner integral is well-defined. 
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Fix t G [0,T], let ^ 1,^2 be two Borel representatives of VPt ,i.e \Fi, Sl /2 are 
Borel functions and T; = VPt. C? — a.e. Concerning the first term, consider the set 

E = {(x,x) G f^ 0 |£('Fi (x))d t if(t,x) ± £(>F 2 {x))dtip(t,x )} 

Then a t (E ) = ah.({x G n o o|’l , i(20) 7 ^ 'I f 2 (al))}) = 0. The argument for other 
terms are the same. □ 

Now we define a more general class of initial data, and state the existence result 
of renormalized solutions with such data. To prove the existence, we need to study 
the functional E v (h , 7 ) and J^(P, R ) with more general v than considered in section 
3. This is done in subsection 2.2. As a byproduct of our proof, we will get the dual 
space existence result with a general measure valued initial data whose support is 
contained in R 2 x [—^, —5] for some S > 0 

Definition 4.4. Let Pq : floo -> R be a convex function, and ho G L 2 (Ll 2 ). Then 
we say ( Po,ho ) is a generalized data if the following holds: 

(i) VP 0 G L 2 (fl ho ) and < d X3 P 0 < -5 

( ii) ho > 0 and / hodx\dx 2 = 1 

(iii) P 0 {x 1 ,X 2 ,h 0 (x 1 ,X 2 )) = \(x\ + x%) \/(xi,x 2 ) G 


Next we show that the notion of relaxed Lagrangian solution defined here is 
consistent with weak Lagrangian solution when the measure at is induced by some 
physical flow map, except possibly the inverse map may not exist. 


Lemma 4.5. Let (P, h, a) be relaxed Lagrangian solution such that P G C([0, T]; iL 1 (fl 00 )) 
with admissible initial data. Suppose there exists a Borel map F : [0, T] x V h0 -t 
Lloo, which is weakly continuous in the sense that 

/ F t (x)ijj(x)dx —> / F to (x)if(x)dx Vif G C&(R 3 ;R 3 ) as t —»• to € [0, X 1 ] 

such that at = (id x F t )^tJh 0 - Then ( P,h,F ) is a weak Lagrangian solution except 
possibly ( iv ) of Definition 3.7. 

Proof. First we observe that the assumption P G L°°([0, T], iF 1 ’ 00 (n oo )) p| C([0, T], iJ 1 (n oo )) 
implies that P G C([0, T\, kF 1,r (f2 00 )) V?’ < 00 . Also observe the following 

/ \F t (x)\ 2 dx = I \x\ 2 dx —» / \x\ 2 dx = / \F to (x)\ 2 dx as t t. 0 

JQhr* J Qh(t ,•) J£lh( + ~ J£lhr. 




The equality above is because 


— projx&t — b;Aj 

and the limit happens is because of the assumption h G C([0, T\ x CI 2 ). Since F t is al¬ 
ways bounded, this combines with weak continuity implies F G C([0, T], L r (crh 0 C 3 )),Vr < 
00 

Now it only remains to check that for Z(t,x) = VPt(F t (x)) (well-defined thanks 
to the above observed F t $<Jh 0 = cr^p,.)), it satisfies (u), but this is seen by taking 

£(y) = Uk in (4.1). □ 


Here is the main result of this section. 
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Theorem 4.6. Let (Po,ho) be generalized data in the sense of Definition 4-4, then 
there exists a relaxed Lagrangian solution (P,h,a) having (Po-ho) as initial data. 
Besides, we have the following continuity in time: 

(i) heC(lO,T};L 2 (D 2 )) 

( ii ) f(P,VP)a h G C([0,T];L 1 (n oo )) G C b (R x R 3 ) 

4.2. Measure-valued solution in dual space. Recall that when v$ e P ac (R 3 ), 
the system in dual space can be written as 

dtv + V • (ww) = 0 

w = J{y -VP*(f,z/)) ye R 3 

(h, (id x VPjjffi) minimizes E v (h,j) 

v\t=o = 

See the end of section 2. 

Since the second equation involves VP* (f, y), it is not well-defined if v is singular,we 
need to find an substitute for VP*. A natural candidate for this is to use barycentric 
projection. 

Definition 4.7. Let /z, v G P(R 3 ),A G T(p,,v), we say the Borel map (defined 
v — a.e a,e ), 7 : R 3 —>• R 3 is the barycentric projection of A to v if the following is 
true 

[ i{y)-l{y)dv{y) = [ [ £(y)-xd\(x,y) G C 6 (R 3 ;R 3 ) 

J R 3 J J R3 X R3 

Now if P t be a family of convex functions and set u t = VPtj(T/ l ( t> .) , then we denote 
7 t(y) to be the barycentric projection of (id x VPt)ja 7 onto 17 , or equivalently 

[ €(y) ■ 7 t(y)dv t (y) = [ £(VP t (x)) • xa h (x)dx G C b (R 3 ; R 3 ) 

It‘s easy to see when v t « £ 3 ,then VP t *(y) = 7t(y) 17 — a.e 

Thus in the general case when 17 is not necessarily absolutely continuous, a 
natural way to write the system in the dual space is (where [ 0 ,T] 9 t 1 —> Z 7 is a 
Borel family of measures) 

(4.8) d t v + V • (nw) =0 in [0, T\ x R 3 

(4.9) W (t, y) = J(y — 7 t(y)) in [0, T] x R 3 

(4.10) (h, (id x V Pt)^cFh) minimizes E v (h,j) 


(4.11) v \ t =0 = ^0 

For immediate use, we first prove a lemma. 

Lemma 4.12. Let v be a solution of above system, fix t G [0,T] and 7 := (id x 
VP b ) jjOTj 6 e f/ie optimal measure with marginals ah. and v t under quadratic cost, let 
7 1 be the barycenter projection defined in Definition 4- 7, t/ien 
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Proof. Observe that 

\\j[id — 7t]||i, 2 ( i/ t ) = sup J(y~it(y))<i>(y)dvt(y) 

0ec c (R 3 ),||0|| i 2<i J 

But 

J J(y - it{y))<t>{y)dv t {y) = J J(VP t (x) - x)^ P t {x))da h ^.){x) 
<(j(xi- d 1 P t (x)) 2 + (x 2 - d 2 Pt(x)) 2 da Ht ,.)(x))?(J 4> 2 (S/P t (x))da h ^.))^ 

f iVPtixtfdx^M]^) 


< [2(max |x|' 
r 2 o 




But since VP^ayo..) = ^t, we know 

[ \VP t (x)\ 2 dx = W 2 (is t ,6 0 ) 


□ 


We denote by AC p (0, T; P 2 (R 3 ))(See also [3] chapter 8 ) to be the set of all paths 
fi : [0,T]3l->/i ( e V -2 (R 3 ) for which there exists /3 G L p { 0, T) such that 

W 2 (ns,lit) < f P{r)dT Vs,t G [0, T] 

J S 

First we want to study the stability of the measure valued dual space solution under 
perturbations of initial data. 

Proposition 4.13. Suppose {uo,uo}n>i C P 2 (R 3 ) with supp is q,is 0 C R 2 x 
[—j, —c5] and such that 

Vq —> iso narrowly and sup M 2 (isq) < oo 

n> 1 

Let is n G AlC' oo (0, T; P 2 (R 3 )) be solutions to the dual space system corresponding 
to initial data isff, and also that supp is™ C R 2 x [— — <5]. Then a subsequence is™ 
converges to a solution v G AC oo (0,T;^(R 3 )) with initial data iso,i- e 

(i) W p {is?,is t ) —> 0 for all t G [0,T], and all p G [1,2) 

(ii) is G AC°°(0,T;V2(H 3 )) is a solution to dual space system with initial data 
no 


Proof. First we have for each n, and 0 < s < t < T 
(4 - 14) 

W 2 (is™,is™) < f || - 7 ”|| L 2 ( ^n ;R 3 )dT < 2 [ W 2 (is™, S 0 )dr + 2 max |x|(t - s) 

Js T Js “2 

The first inequality used [3] Theorem 8.3.1. The second inequality is by Lemma 
4.12. 

Now take s = 0, we then have 

W 2 (< 5 0 , K) < sup W 2 (6 0 , is™) + W 2 (i/J, is™) 

n> 1 


< sup W 2 (So,isq) + 2 

n> 1 


W 2 (is™, §o)dr + 2max|a;|t 
n 2 


o 
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Then we can apply Gronwall to obtain 

(4.15) W2(5q, vf) < siip W 2 (Sq, Vq )e 2T + (e 2T — l)Tma x |a;| 0 <t<T 

n> 1 

which means {vt}te[ o,T],n>i are uniformly bounded in T-^R 3 ), so it is tight. 
Combining (4.14),(4.15), we get 

(4.16) W 2 (t'",i'") < 2 e 2T [sup W 2 (So,Vq ) + (T + l)max|x|]|t - s| 

n> 1 n 2 

which means {nfj n >i are equi-continuous in t under the distance W 2 . This allows 
us to take a subsequence (not relabeled), such that 

v " —> v t narrowly Vt G [0,T] 

Since sup n>1 f \y\ 2 dvf < 00 , we see 


J \y\ p dv? ->• J \y\ p dvt Vp G [ 1,2) 


this implies W p (v^,v t ) —> 0 by [3] Proposition 7.1.5. Also it‘s easy to see v t G 
AC°°(p,T;V 2 ( R 3 )) by (4.12). Finally let's check v t solve the dual space system. 
Take <j>(y,t) = x(t)v(y) be test function such that x{T) = 0, then for each n, we 
have 


/r 3 


x'{t)v{y)dvt{y) 


/ R 3 


J(y- 7 t(y)) ■ Vii(y)dv?{y) 


+ [ x(Q)v(y)dv£{y) = 0 

J R 3 


where 7 f(y) is the barycentric projection of (id x VP")^/,™ onto 1 /". 

Because of narrow convergence, and W v convergence for p G [1,2) already 
noted,we see / Q T J R3 x'(t)ri(y)dv?(y),f R3 x(0)r/(y)distf (y) ,/ Q T f R3 Jy ■ S7y(y)duf(y) 
will converge to the right limit,while 


/ 0 Jr 3 


Ji t(v) ■ Vr](y)dis?(y) = 



Jx • V77(VP"(a;))crftn (x)dx 


Notice that above integrand does not involve £ 3 , by the definition of J, and Xi,X 2 
are bounded, we can conclude by using Theorem 2.53(m) that 



Jx ■ S7r](yPt(x))(Jh™(x)dx —> 



Jx ■ Vp(V Pt(x))<Jh t {x)dx 


= [ [ Jjt(y) • Vy(y)dv?(y) 

Jo J R 3 

□ 


As a corollary to above proposition, we can deduce the existence of measure¬ 
valued solution in dual space with possibly unbounded support. 

Corollary 4.17. Let i/q G ^(R 3 ) and such that supp uq C R 2 x [—^,—<5], then 

there exists a family of measures [0, T] B 1 1 —> vt such that v G AC°°( 0, T ; ^(R 3 )),supp u t C 

R 2 x [—-y, —<S] and solves the dual space system. 
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Proof. Define uj = „ o[B l 0 ^ n)) WoXb( o,n)] *j±, then v o e ^(Rs), supp C R 2 x 
[—4, —5] and compact. Also we have Ug —> narrowly. Then we can use theorem 
2.30 to get a weak solution uf in dual space solution. The argument used to show 
(4.12) shows u” e AC oo (0,T;P2(R 3 )) and has Uq as initial data with supp v ” C 
R 2 x —5]. The previous proposition gives us a solution v t with initial data vq 

□ 

4.3. Existence of relaxed Lagrangian solution with generalized data. In 

this section, we will prove the existence of relaxed solutions with generalized initial 
data. First we show that for the generalized data (Po,ho) defined previously, Poln^ 
gives a maximizer. This is a generalization of Proposition 3.3. More precisely, we 
have 


Lemma 4.18. Let ( Po,ho ) be a generalized data, put vq = VPojjCpg,. Let (Po,Po) 
be a pair of maximizer of J Uo (P, R), which satisfies (3.39), (3.)0) , then Pq = Po on 
Llh 0 and (ho, (id x VPo)pcr^ 0 ) is the minimizer of E^ 0 (h,^) 

Proof. The argument is similar to Proposition 3.3. The difference is that one need 
to define 

P 0 (x) = max(Po (x) , ~(xj + x*)) 

and 

R 0 (y) = sup (x-y- P 0 (x)) 

XdzQoo 

This definition ensures Po is finite. The rest of the argument works in the same 
way as Proposition 3.3. □ 

Remark 4.19. Recall Definition 4.2, the same argument as in Proposition 3.3 shows 
that for each fixed t £ [0, T], if we put v t = V Pt^hft,-)-, then (h t , (id x VP t )$ah t ) 
is the minimizer of E Vt (h,^), Pt is the restriction of some maximizer of J Vt (P,R) 
restated on f \ t 

First we observe that, similar to the fixed boundary case, the relaxed Lagrangian 
solution will give rise to a measure valued solution in dual space. 

Proposition 4.20. Let (P, h , a) be a relaxed Lagrangian solution and let vt = 
VPtijOT,^.), then 

(i)v solves the dual space system with initial data vo = VPojaAg, 
(%e4C”(0,r;P 2 (R 3 )) 


Proof. We choose test function as <j>(t,y) = x(t)^(j/), then we can compute 


/ o J R 3 


X'(t)il>(y)dv t (y)dt = 


X' (t)i>('VP t (x))da h ^ ! .)(x)dt = 


i’(VP t (x))x'(t)da t (x,x) 


= -/ / VV'(VP t (i)) • J(VP t (x)-x)x(t)da t (x,x)dt- / ip(VP o (x))x(0)dx 

Jo Jn 2 ^ Jn hQ 

In the above, we have used (4.1) and the marginal properties of dat(x,x). Now 


Vf’(VP t (x))-J(X7P t (x))x(t)da t (x, x) = 


Vij}(VPt(x))-J(VPt(x))x(t)dxdt 




V^(y) • Jyx{t)dv t {y)dt 
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while the term 

r T r r T 

\ 7 ip(VP t (x))-Jxx(t)dat(x,x)dt = 


/0 J f2L 


/o 


\/'ip(VPt(x))-Jxx(t)da Htt .)(x)dt 


The last term 


= / / W(») • J^ t {y)x{t)dvt{y)dt 

Jo JR 3 

V>(VP 0 (a;))x(0)ote = [ tp(y)x(0)du 0 (y) 




/ R 3 


So we proved (i). To see (n) is true, we notice 


v T )dr < / 2(max|a;| + do))cP 

./« ^2 

In the above, we used Lemma 4.12. Now we can use Gronwall inequality to conclude 
that {z'tjteio.T] is bounded in P 2 (R 3 )- □ 

The existence result of relaxed solution is an easy consequence of the following 
lemma. 


W 2 {v s ,v t ) < J \\J(id — 7 t )|| l 2 ( 


Lemma 4.21. Let (P 0 n , ),(Po, ho) be generalized initial data in the sense of Def¬ 

inition 4-4> Suppose also that 

sup||VP 0 n ||i 2 (n ) < oo 

n> 1 0 

and 

K -»• ho in L l (D 2 ) Z{Pg,VPg)<r h » -> C(Po, VP 0 K in L 1 ^) 

VC € C b (R x R 3 ) 

Let ( P n ,h n ,a n ) 6 e relaxed Lagrangian solution with initial data {Pff,h q), f/ien 
t/iere exists a subsequence which converges to a relaxed solution with initial data 
(Po, h 0 ). More precisely, for any t G [0, T] 

( i ) h " - 4 - h t in L r (Q 2 ) Vr e [1,2); 

(ii) a n —>■ a narrowly; 

(Hi) C(R”, Vi?)®*- -)• C(Pt, VP*K m PHfloo) VC G Cb(R x R 3 ) 

Proof. Define z/” = VP^cr/q*, then by Proposition 4.20, we know z/ n solves the dual 
space system with z/ n G AC' oo ( 0 ,T;P 2 (R 3 ))- Since —| < c^Pt < —5, we also 
know that supp vf C R 2 x [—|, — 5] 

Define Z 2 0 = V Po;<Xh 0 . By assumption, one has supp v r f ■ ,y o C R 2 x 
Also 

[ y(y)dvo (y) = [ r]('VP^(x))a h n(x)dx -> [ r/(V P 0 (x))a ho (x)dx 
Jr 3 Jn oc dOoo 

= [ y(y)dv 0 (y) e c&(R 3 ) 

•/R 3 

So z 2 q —» z^o narrowly. Also we have sup ra>1 M 2 (vq ) < oo, because we assumed 
sup„>i ||VP 0 "||i 2 (nji „) < oo. 

We are now in a position to apply Proposition 4.13 to get a subsequence of vf 
which converges to a solution v G AC°°(0, T; P 2 (R 3 )) which solves the dual space 
system with initial data vq. 

By (4.11) in the proof of Proposition 4.13, we know {^"}zg[o,T],n>i is uniformly 
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bounded in V 2 (R 3 ), hence {||VP t "|| L 2 (n K) } tG[0 ,T], n >i are uniformly bounded. Re¬ 
mark 4.19 shows for each t G [0, T],(/i", (id x VP")}^™) is the unique minimizer of 
E v n(h, 7 ). Now let (P", P"),(P, P) be maximizers of J I/ n(P,R),J Vt (P,R) respec¬ 
tively satisfying (2.39),(2.40), then by Theorem 2.53, we conclude for each t G [0, T\ 

h™ —> h t in L r (Q 2 ) for any r G [1,2) 

£(P t ", VP t ")a fc « -»• £(P t ., VP t ) T 1 (fi c >o) V? G C„(R x R 3 ) 

Also we know that {||/i"||L 2 (n 2 )}ie[o,T],ri > 1 is bounded by a universal constant 
by Lemma 2.46. By Remark 4.19, we actually have P” = P” on O b n, so we really 
have 

£(P t ", VP”H ? -A £(P t , VP,K in L\n oo) G C b (R x R 3 ) 

Finally since for each n, proj x a " = ah^,pr°jx a t = a K- Because of the L 2 bound 
of hf, we know {a 11 }, as a measure on [0, P] x x floo, is tight. Therefore up to a 
subsequence, we can assume a n —> a narrowly, and the limit a must disintegrate as 
da(x,x) = dat.dt , with proj x at = <Jh 0 and projxQ-t = 0 h t - So far we have checked 
point (i)-(iii) in Definition 4.2. (v) follows from the assumed convergence of the 
initial data and the fact that Pff —> P 0 locally uniformly on x {0}. It only remains 
to check (P, h, a) satisfies point (iv). For each fixed n, we have 




^CVP[ l (x))dttp(t, x)daf(x , x)dt+ 


r 

0 



V^(VP t n (£))-J(VP t n (i)-i)^(t, x)do$(x , x)dt 




42 


JINGRUI CHENG 


First we estimate B. We notice that 

[ [ (£(VP t (x)) - ^{g(t,x)))d t ip(t,x)da n (t,x,x) 

Jo 

(4.24) < ||^||l» [ T [ \£{VP t (x)) - £(g{t,x))\a h n(x) 

Jo ifioo 

<\\dtJ>\\L°°[ [ f \€(VPt(x))-t(g{t,x))\a h »(x)+\\Z\\ L 'fJ f (h?-K)+d Xl dx 2 ] 
JO J^lur JO J Clo 


< IMl~[||£||l ~6 + ||^|| ioo 11^(02)] < 2||^|| i oo||£|| L «e 

The same estimate also works for D. while C —> 0 by narrow convergence. 
It only remains to estimate A. 

[ [ (£(VP"(x)) -f;(\7Pt{x)))dtip(t,x)da n (t,x,x) < 

Jo Jn%, 


WdtipWLoo / |£(VP”(x)) -€{VPt(x))\a h *(x)dxdt 

Jo Jn oo 

Then we can write 

|£(VP/*(x)) -£(VP t (x))K ? (x) < UVP^{3A,)a K {x)-^P t {x))<j ht {x)\+ 

\\t\\L°°Wh?{x) ~ cr ht {x)\ 

By the convergence already noted, their integral goes to zero. So we can pass the 
first term of (4.22) to limit. 

For the second term, the term f Q f Q2 V£(VP"(x)) • Jxip(t, x)da^(x, x)dt can be 
passed to limt in the same way as first term, if one notices that Jx does not involve 
£ 3 . We only need to deal with f Q2 V£(VP t n (x)) ■ JVP t n (x)if>(t,x)dat(x,x)dt 

Let H e (^(R 3 ) be a cut-off function such that H(y) = 1 if |y| < K, H(y) = 0 if 
| y | > K + 1, K to be determined below. Then we can write 

(4 t 5) 

f f VZ(VP?(x))-JVP?(x)dtil>dc%(x,x)- f I V^{yP t {x))-JVPt{x)dt^da t {x,x) 

Jo Jn? Jo 


<II^IU«||veilL- 


|VP”(x)(l - H(VP t n (x))\a h n(x)dx 


to 


to -tn o 


|VP t (x)(l - H{VP t (x))\a ht (x)dx] 


+ [ 


to Jn’L 


V£(VP t n (x)) • JVP?(x)H(VP?(x))dtil>do£{x,x) 


- f I V£(VP t (x)) ■ JVP t (x)H(VP t (x))dtHat(x,x)\ :=A + B + C 
Jo J 

For A , since c” = VP^jjcr^, we can choose K large, so that 


A = 


to J R 3 


I VO ~ H(y)\dv?(y) < 


to J{\y\>K} 


I U n A , s Tsn Pn,t M 2(^t) 

\y\dv™dt <- —7 - < e 


A' 
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The term B is exactly the same as above. For term C, it can be dealt with in the 
same way as first term of (4.22) already shown above, since the integrands in C is 
bounded because of the cut-off. □ 

Now we are ready to prove Theorem 4.6, the existence of relaxed Lagrangian 
solutions. 


Proof. Let v 0 = VPojoTk,, then u 0 G ^(R 3 ) and supp n 0 C R 2 x [— y, — <5]. Define 
u 0 ” = Vo tB(o,n)) [Xb(o, 71 )^ 0 ] *j a- Then z/q has compact support and smooth and 
Up —> uq narrowly. Let (Pq,Rq) be the maximizer of J^n n (P 1 R) with H n taken 
large enough. (Kq, 7 J) be the minimizer of E v n(h, 7 ), then (Pq%/iq) are admissible 
initial data. We know by Lemma 2.54 that it is also the maximizer of J„»(P, P) 
after suitable extension of Pq, hence they are also generalized data in the sense 
of Definition 4.2. Also by Theorem 2.53, we have fift —>■ ho in L r (Pl 2 ), for any 
r G [1,2), and £(P",Vi^H» £(P 0 ,VP 0 W 0 in L 1 ^). Let (P n ,h n ,a n ) be 

the Lagrangian solution with initial data ( Po,hfj ) given by Theorem 3.9. Then 
previous lemma gives us a relaxed solution with initial data (Po,/io) 

The continuity property in time of h is given by Theorem 2.53(iv) upon noticing 
that v G AC°°(0,T; P 2 (R 3 )) ' □ 


5. Appendix 

Here we prove that 

Lemma 5.1. Given (Po, Rq) be a pair of convex conjugate maximizer of Jff (P, R), 
let ho = hp o be such that 0 < ho < H, Suppose also Q 2 C Bd(0), where A = 
B d { 0) x define 

Pi (y) = sup [x • y - max(P 0 (x), \{x\ + x%))\ 

Pi(x) = sup(x-y-Pi(y)) 
v &a 

Then the following holds 

(i) (Pi,Pi) are convex conjugate over SIh, A 

(ii) Pi(xi,X2,0) > |(x 2 + x%) V(xi,X 2 ) G f h 
(m)Pi(xi, X 2 , 0) < ^(x 2 + x$) whenever ho(x±,X 2 ) = 0 
(iu)(Pi,Pi) is a maximizer of J^f (P, R) and /ij^ = ho 

Proof. First we prove (i), denote Po(x) = max(P 0 (x), \{x\ +x\)), we need to show 

Ri{y) = sup (x • y Pi(x)) 

x€lQ,h 

By definition of Ri(y ), we know 

Ri(y) + Po(x) > x • y Vx G SIh , Vy G A 

Hence 

p 0 (x) > Pi(x) x G pin 

Ri(y) < sup (x-y-Pi (x)) 
x£Q,h 


So 
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By the definition of Pi, we have 

Pi (a:) + Ri(y) > x ■ y Vx G f Ih ,Vy G A 
So 

Ri{y) > sup (x-y-Pi(x)) 
x££Ih 

(i) is proved. Now we prove (ii). We can observe that 

Ri{y) < sup (x-y-\(x\ + x 2 2 )) 
x£Qh ^ 

Fix (x5, x°) G by assumption, we can find zq G [—4, —<5], such that (xj, x°, 20 ) G 

A, then 

Rl(Xi,X 2 ,Z Q ) < sup (X- (Xi,X 2 , 2 0 ) - ^(*1 + X 2 )) 
x£Qh " 

= sup(xiX? + x 2 X2 - \{x{ +X2)) < ^[(£C?) 2 + (a^'a) 2 ] 

o 2 2 2 

In the equality above, we noticed X 3 Z 0 < 0. Therefore, 

Pl(xi, x°, 0) = sup(x?yi+x^ 2 -Pi(y)) > (x?) 2 +(x£) 2 -Pi(x?,x£, z 0 ) > ^[(x?) 2 +(x£) 2 ] 

(ii) is proved. Now we prove (Hi). Fix (x^Xg) G such that ho(x\,x®) = 0, then 
(x?,x§,o) ^ n ho 

But 

Pi(y) = max[ sup (x ■ y - P 0 (x)), sup (x • y - P 0 (x))] 
xe£ih 0 x£Q. ho 

On tt ho ,P 0 (x) > \(x\ + Xj), so P 0 (x) = P 0 (x), otherwise P 0 (x) = ^(xf + x|) 

Hence 


Then 


Pi(i/) > sup (x • y - -(x? + x^)) 

X^hr, ^ 


(x°,x°, 0) ■ y — Ri(y) < -((x?) 2 + (x°) 2 ) y G A 


Taking supremum over y, we obtain 

Pl(x?, x°, 0) < i((x?) 2 + (x°) 2 ) 

(Hi) is proved. Finally we prove (iv). We start by showing that Pi = Po on ilf l0 
Indeed 


Ri(y) = sup (x-y- P 0 (x)) < sup (x • y - P 0 (x)) = P 0 (y) 
x^Qh x£Q,h 


So 


Pl(x) = sup(x • y - Pi(y)) > sup(x • y - Ro(y)) = Po(x) 

yCA yCA 

On the other hand 

Ri(y) > sup (x-y- P 0 (x)) 

x£flh 0 


Pi(x) = sup(x • y - Ri(y)) < Po(x) 

y€A 


So if x G H/j, 
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Combined with ( ii ), (Hi) above and recall the definition of hp , we can deduce that 
hp o = hp x = ho- Hence if we define 70 = (id x VP 0 )jjO7 l0 , then (ho, 70) is the 
minimizer of E v (h, 7), So 


J?(Pi,Ri) 


j [\{vl + yl) - Ri{y)\dv(y) 



— P\(x)]dx 


s 

J 1 


[^( 2 /i +vl) + \{x\+x 2 2 ) - Ri(y) 

xA z z 


Pi(a;)]d 7 o = E„(h 0 , 70) 


Since (Pi, Pi) convex conjugate, and Po = Pl on flh. 0 , R\(y ) + Pi (x) = x ■ y 70 — 
a.e □ 


Remark 5 . 2 . We notice that above argument still works even if H = 00 


Lemma 5.3. Suppose (Pi,Pi) are maximizers of Jff (P,R), convex conjugate over 
VLh and A. Suppose also that Pi(xi,x 2 , 0 ) = ^(xf + x 2 ) whenever hi(xi,x 2 ) = 0 , 
where hi = hp < H. We define 

R2(y)= sup (x-y-Pi(x)) 

x€fl hl U{^3=0} 


P 2 (x) = sup(a: • y - R 2 (y)) 
yeA 

then 

(i) Pi = P 2 on fl hl U{a.’3 = 0 } 

(ii) (P 2 l R 2 ) are convex conjugate over fiu and A. 

(iii) (P 2 , R 2 ) is also a maximizers and hp o = hi. 


Proof. Take x £ f 2 /^ IJ{- r 3 = 0 }, then \/y £ A, we have 


x ■ y - R 2 (y) = x ■ y - sup (x ■ y - Pi(x)) < Pi(x) 

U{®3=0} 


Take supremum over y to get 

Pi(x) < P\ (x) (Vcc € Ct hl l^J{cc 3 = 0}) 


On the other hand, by definition R 2 (y) < Ri(y), since (Pi,Pi) conjugate over fin 
and A, we obtain 

Pi (a;) < P 2 (x) (Vx £ Qh) 

(i) is proved. 

To see (P 2 ,R 2 ) also convex conjuate over f 1 h and A, we only need to show 


R 2 (y) = sup (x-y- P 2 (x)) 
xEQh 

Obviously LHS < RHS by (i). 

On the other hand, for all x £ Qh 


x-y- P 2 (x) =x-y- sup (a; • y - R 2 (y)) < R 2 (y) 

yeA 

So (ii) is proved. 

Finally we only need to see jff(Pi,Ri) = Jff (P 2 , R 2 ). By assumption, we know 
(hi, (id x VPi) # <r fcl ) is the minimizer of E v (h, 7) 

By (i), we know that 

P 2 (xi , x 2 , 0 ) = — (ajf + x|) whenever hi(xi,x 2 ) = 0 
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Combining the fact that Pi = P 2 on f lh x , it‘s easy to see hp 2 = hi. The same 
argument as previous lemma (iv) shows that (P 2 ,i? 2 ) is a maximizer. □ 

We derive the following corollary as an easy consequence of previous two lemmas. 

Corollary 5.4. Let A = B d(0) x <5] be such that fl 2 C B d(0) ; and H is 

chosen such that J^f(P,R) has convex conjugate maximizers over Tin and A, say 
(Po,Ro), and ho := hp o < H. Then there exists a maximizer (P 2 ,i? 2 ) of jj?(P, R) 
which satisfies the following conditions: 

(i)(P 2 ,P 2 ) are convex conjugate over both Llh 0 = 0}, A and LIh,A. 
(u)P 2 (;ri, a; 2 , 0) = |(a;f + x%) whenever ho(xi,X 2 ) = 0. 

Proof. Let (Pi, Pi) be the pair given by Lemma 5.1. The conclusion of Lemma 5.1 
shows (Pi,Pi) satisfies the assumptions of Lemma 5.2. Let (P 2 ,P 2 ) be the pair 
given by Lemma 5.2. Then such a pair is a maximizer by Lemma 5.2 (Hi). They 
are convex conjugate over Lth 0 = 0} by their very definition. They are convex 

conjugate over Qh, A by Lemma 5.2 ( ii ). P 2 satisfies (ii) because of Lemma 5.2(*) 
and Lemma 5.1 (ii), (in). □ 

The following lemma can be found in m section 5.3 Theorem 1, and so we omit 
the proof. 

Lemma 5.5. Let fl be a convex domain in R , let P : fl —> R be a convex function 
such that ||P||_Li(n) < 00 . Let fli = {x £ fl\dist(x, fl c ) > r}. Then there exists a 
constant C = C(\\P\\i,i{Q.),r,d),such that 

(5-6) ||P||L~ ( n 1 )<C' ||VP|| L ~ (ni) <C 
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